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NQnlinear System^ Method o£ Design thereof andL — 
Corogulier Program ggod-uct 

The present: invention relates tic non-linedr systems, 
mechods of design thereof in the frequency domain and 
computer program products. More pa.rticulariy^ the present 
invention relates to a non-linear system having pre- 
determlnable frequency response characteristics- The 
invention can be utilised to design and realise, for 
example, nonlinear 'filters having a required frequency 
response or transfer functions having specified transfer 
characteristics or within a control system context. 



The possible frequency components in an output signal 
of a linear sysitem are exactly the saitie as the frequency 
components of a corresponding input signal. Conventional 
linear filter design is based on the principle that energy 
in unwanted frequency bands is attenuated. 

The Dolby filter^ which varies the amplitude of the 
output signal as a function of the level and frequency of 
the inputs is an example of a nonlinear filter system. 
However/ when compared with the inputs the output does not 
contain any additional f recjuency components - Modulation is 
another concept related to nonlinear filtering/ which is 
associated with signal transmission where the signal to be 
transmitted is modulated by a carrier signal and then 
transmitted through a medium. Although a modulation device 
allows energy, to be moved from one frequency band to 
another, the output frequency components of such a device 
depend not only on the input components but mainly on the 
carrier signal. Therefore, the energy transfer implemented 
by modulation is realised by a two input and one output 
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system where one inpur is rhe carrier signal and the other 
input is tne signal ro be processed. 

The prior art lacks a non-linear sys-tem ^nd method 
5 of /apparatus for the design of such a non-linear system for 
predictably transferring energy from one freciuency or band 
of frequencies of an input signal to another frequency or 
band of frequencies independently of any other input signal* 
Further, the prior art lacks a non-linear control system 
10 which can predictably transfer energy from one frequency 
band to another frequency band. 

It is an object of the present invention to at least 
mitigate some of the problems of the prior art. 

15 

Accordingly, a first aspect of the present invention 
provides a method for -designing a non-linear system for 
transferring energy from a time or spatial domain input 
signal having a first spectruiri at a first pre-dererminable 
20 frequency or range of frequencies to a time or spatial 
domain output signal having a second spectrum at a second 
pre-determinable frequency or range of frequencies. 

Preferably, the method comprises the steps of 

identifying the first spectrum of the time or spatial domain 
2S input signal from which energy is to be trans ferred, 

specifying the second spectrum of the time ox spatial domain 
output signal to which said energy is to be transferred; and 

calculating, using a frequency domain description of said 
output signal/ for example, the output spectrum^ expressed 
30 in terms of a frequency domain description of said input 
signal and coefficients of ' a time or spatial domain 
description of a generalised uon-linear syst^m^ the 
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coefficients of a said time or spatial domain description of 
said generalised non-^linear system in order to give effect 
to the energy trans fer- 



allQws energy at -a particular frequency within a given 
system to be transferred to another frequency or band of 
frequencies at which the response of the system is greatly 
reduced or negligible,- or 

10 allows energy of a signal which is transmitted at a 
particular frequency or band ' of frequencies to be 
transferred, without using an additional modulating signal, 
to another frequency or band of frequencies at which the 
associated transmission media allows signals to pass, or 

15 allows energy at a particular band of frequencies to be 
transferred and spread over a new wider range of frequencies 
so as to attsiaiiate the energy by employing the desired 
interkernel and intrakernel effects of nonlinear systems* 

20 The present invention is based upon the relationship 

between the input and output spectra or freqfuency components 
of nonlinear systems / and the relationship between the input 
and output frequencies and/or frequency ranges in the 
nonlinear case. In addition, the invention utilises a 

25 mapping between the time or spatial domains and frequency 
domain which allows the ourput spectra or frequency content 
of nonlinear systems t:o be described completely by- the 
coefficients of time or spatial domain models which 
represent the filter or non-linear system to be constructed, 

30 A second aspect of the present invention provides a 

method for manufacturing a non-linear system for 
transferring energy from a time or spatial domain input 



5 



Advantageously^ the present invention 
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Signal having a first spectrum at a first pra-detaTnninable 
frequency or range of frequencies to a time or spatial 
domain output signal having a second spectrum at a second 
pre-determinable frequency or range of frequencies, said 
method ot manufacture comprising the steps of 

(a) designing said non-linear system comprising the steps of 
identifying the first spectrum of the time or spatial domain 
input signal from which energy is to be transferred, 
specifying the second -spectrum of the time or spatial domain 
output signal to which said energy is to be transferred, and 
calculating, using a frequency domain description of said 
output signal, for example, the output spectrum, expressed 
in terms of a frequency domain description of said input 
signal and coefficients of a time or spatial domain 

15 description of a generalised non-linear system, the 
coefficients of a time or spatial domain description of said 
generalised non-linear system in order to give effect to the 
energy transfer, and 

(b) materially producing the non-linear system so designed. 



10 



20 



A third aspect of the present invention provides a data 
processing system which can transfer energy from a time or 
spatial domain input signal having a first spectrum at a 
first pre-determinable frequency or range of frequencies to 
25 a time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or range of 
frequencies, said system comprising 

means for identifying the first spectrum of the time or 
spatial domain input signal from whxch energy is to be 
30 transferred/ 
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means for specifying the second spectrum of the time or 
spatial domain output signal to which said energy is to be 
transferred/ and 

means for calculating, using a frequency domain description 
5 of said output signal, for example, the output spectrum, 
expressed in terms of a frequency domain description of said 
input signal and coefficients of a time or sparial domain 
description of a generalised non-linear system, the 
coefficients of a tictie or spatial domain description of said 
10 generalised non-linear' system in order to give effect to the 
energy transfer, 

A fourth aspect of the present invenrion provides a 
computer program product for designing a non-linear system 

15 for transferring energy from a time or $patial domain input 
signal having a first spectrum at a first pre-determinable 
frequency or range of frequencies to a time or spatial 
domain outpu-t signal having a second spectrum at * secjorid 
pre-decerminable frequency or range of frequencies, said 

20 computer program product comprising 

computer program code means for identifying rhe first 
spectrum of the time or spatial domain input signal from 
which energy is to be transferred, 

computer program code means for specifying the second 
25 spectrum of the time or spatial domain output signal to 
which said energy is to be transferred, and 

computer program code means for calculating, using a 
frequency domain description of said output signal, for 
example, the Qutput spectrum, expressed. in terms of a 
30 frequency domain description of said input signal and 
coefficients of a time or spatial domain description of a 
generalised non-linear system, the coefficients of a time or 
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spatial domain desc^ription of said generalised non-linear 
system in order to give effect to the energy transfer. 

A fifth aspect of the present invention provides a noa-* 
5 linear system which can transfer energy from a time or 
spatial domain input signal having a first spectrum at a 
first pre-determinable frequency or range of frequencies to 
a time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or range of 
10 frequencies/ said system comprising 

means for identifying the first spectrum of the time or 
spatial domain input signal from which energy is to be 
tirana f erred, 

];p_ea_jXS for specifying tiie second spectrum of the tim^ or 
15 spatial domain output signal to which said energy is to be 
transferred, and 

means for giving effect to the energy transfer using 
coefficients of a time or spatial domain description of a 
generalised non-linear system, said coefficients having been 
20 calculated using a frequency domain description of said 
output signal, for example, the output spectrum, expressed 
in terms of a frequency domain description of said Inpxit 
signal and coefficients of a time or spatial domain 
description of a generalised non-linear system, 

25 

Ad-sraiitageously, the fifth einbodlraent allows processing 
£Qr the determination of the coefficients to be performed 
off-line and merely incorporated into a non-linear system 
which uses the coefficients. 

30 
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Embodiments of the presenr invention will be described 
by way of exainples only, with reference to the accompanying 
in which: 

figure 1 shows the effect of traditional signal 
5 processing by^^ for example, a linear filters- 
figure 2 illustrates signal processing according to one 
aspect of the present invention; 

figure 3 depicts a further example of the signal 
processing according to the present invention; 

10 figure 4 shows a non-linear system arranged to give 

effect to the energy transformation shown in figure 3; 

figure 5 illustrates a further energy transformation in 
which energy is distributed over a wider frequency band; 

figure 6 illustrates the power spectral densities for 
15 the input and output signals of a designed nonlinear system.; 

figure 7 illustrates the power spectral densities for 
the inpur and output signals of another designed nonlinear 
systems- 
figure 8 depicts schematically a non--linear system; 

20 figure 9 shows a digital implementation of the non- 

linear sysnem shown in figure 8; 

figure 10 illixstrates the power spectzral densities for 
rhe input and output signals of the nonlinear system which 
was obtained using a design which improves the filtering 

25 effect shown in figure 6; 

figure 11 illustrates the power spectral densities for 
the input and output signals of the nonlinear system which 
was obtained using a design which improves the filtering 
effect: shown in figure 7; 
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figure 12 shows the time doiuain input and output of the 
nonlinear system with the fxrequency domain filtering effect 
shown in figure 10; 

figure 13 shows the time domain input and output of the 
nonlinear system with the frequency domain filtering effect 
shown in figure 11; 

figure 14 depicts the structure of a designed nonlinear 
system; 

figure 15 illus^trates the frequency spectrum of a 
signal to be processed using the present invention; 

figure 16 illustrates the result of a Fast Fourier 
Transform of an input signal having the spectrum shown in 
figure IB; 

figure 17 illustrates the results of the n-dimensional 
(n= 2 and 3) convoluticn integration for the spectxruiu shown 
in figure 16 ; 

figure 13 illustrates the c^utput magnitude frecjuency 
response of a designed nonlinear system to an input signal 
having the frequency spectrum shown in figure 15; 

figure 19 shows the frequency spectrum of a further 
signal to be processed using the present invention; 

figure 20 illustrates the output magnitude frequency 
response of the same nonlinear system as shown in figure 18 
to the further signal to be processed having the frequency 
25 spectrum shown in figure 19; 

figure 21 illustrates the frequency spectrum of a still 
further signal to be processed using the present inventions- 
figure 22 illustrates the output magnitude frequency 
response of the same nonlinear system as in figure IS to the 
still fuirther signal to be processed having the fj^equency 
spectrxim shown in figure 21; 



: 15 



20 



30 
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figure 23 illustrates "ctie structure of another designed 
nonlinear system; 

figure 24 shows the continuous time realisation of the 
discrete rliue system in figure 23; 

5 figure 25 shows at mechanical impleraentation of the 

continuous time system in figure 24; 

figure 26 illustrates the result of a Fast Fourier 
Transform of the further signal shown in figure 1^; 

figure 27 Illustrates the output magnitude frequency 
10 response of a further designed nonlinear system to the 
further signal shown in figure 19; 



figure 28 illustrates the output magnitude frequency 
response of a still further designed nonlinear system to the 
IS signal shown in figure 15; 

figure 29 illustrates the output magnitude frequency 
response of the nonlineair system shown in figure 28 to the 
further signal shown in figure 19; 

20 figure 30 illustrates a flow chart for designing a 

nonlinear system according to an embodiments- 
figure 31 illustrates a flow chart for designing a 
nonlinear system according to a further embodiment; 

25 figure 32 shows the structure of a nonlinear filter 

designed based on specifications for both the magnitude 
and phase of output frequency responses; 

figure 33 depicts the input and output magnitude 
30 frequency characteristics of a specific nonlinear filter 
designed based on specifications for both the magnitude 
and phase; 



9 



SUBSTITUTE SHEET (RUX^E 26) 



30-08-00 11:54 FrorfcHARRISON GODDARD FOOTE 



+441142730312 



T-663 P.ir/99 Job-687 



WO 99/45644 t»CT/GB99/0OSS0 

figure 34 shows the phase angle of the spectrixm 
^zU^) figure 32 in the specific design case shown in 
figure 33 which reflects the phase response 
characteristic determined by the designs- 
figure 35 shows the phase cha.racteristic5 of "Che 
linear phase FIR filter in the specific design ca.se shown 
in figure 33/ 



10 figure 3 6 depletes the discrete time model of a 

nonlinear filter designed to focus energy from two 
different frequency bands into a Single frequency band; 

figure 37 shows the spectrum of an input signal o£ 
ffiis the nonlinear filter in figure 36; 

figure 38 shows the frequency response of the 
nonlinear filter in figure 36 to the input in figure 37, 
which indicates an energy focus effect of the nonlinear 
:^20 filter; 

n= figure 39 illustrates the block diagram of a spatial 

domain nonlinear filter; 

25 figure 40 depicts the power spectral densities a£ an 

input and the corresponding output of the spatial domain 
nonlinear filter in figure 39; 



figure 41 depicts an spatial domain input and 
30 corresponding output of the filter in figure 39; 

figure 42 shows a one dimensional image to be 
processed by the spatial domain nonlinear filter in 
figure 39; 

35 
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figure 43 shows the one dimensional image obtained 
by processing the image in figure 42 using the spatial 
domain nonlinear filter in figure 3 9; and 

5 figure 4 4 illustrates a daca proccsssing system upon 

which erobodiments of the present inTrention can be 
implemented. 



Referring to figure X, there is sH^wn nhe principle otf 
10 for example^ traditional low pass, high pass, and band p^ss 
filtering. Figure 1 shows the power spectrum of a signal 
100 both before and after processing. The energy of a 
signal 100 to be filtered comprises two partSy namely, a 
first part 102 for further processing or of interest and a 
15 second pa.j::t 104 which is of no incerest . Typically^ the 
second part 104 of the signal is attenuated wh-ioh results in 
a second signal 106. The second signal 106 comprises the 
original or a copy of the first part 102 and an attenuated 
portion or an attenuated version of the second part 108, 

20 

Figure 2 illustrates the principle of signal processing 
according to one aspect of the present invention > Figure 2 
shows the power spectruici of a signal 200 both before and 
after processing. The signal comprises a first portion 202 

2S and a second por-tion 204- The first portion 202 of the 
signal 200 is of interest for further processing or output. 
Accordingly^ as a result of the signal processing using the 
present invention, the first portion 202 is retained and the 
energy in the second portion 204 is translated to another 

30 frequency band 20 6. 
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I, DETAILED DECRIPTION 

The theory and method underlying' trtie present invent3.on 

will UQW be described in general terms in steps (i) to (vi) , 



5 (i) Determine the frequency spectruia of a. signal to be 
processed, including the rancte o£ frequencies of the signal. 



10 



15 



20 



(ix) Specify the frequency spectrum of the output signal. 

(iii) Detenuiine the structure of a Noniinea^r Auto-Regressivs 
model with exogenous inputs (NARX model) to ensure that the 
energy transformation between different frequency bands and 
other design requirements, for example, specifications for 
magnitude and/or phase of the output spectrum over the 
required output frequency band can be met or realxsed. 

The general expression for a NARX model is given by 



N 



where y„{k> is a *NARX nth-order output' given by 

(2) 

with 

25 p + q = n, 1- =1,.,.,K^, i = l,,.,,p + q, and X - S'^' X 

K„ is the maximum lag and y(.). and c^^{ -) are the 

output^ input, and model coefficients respectively, A 
specific instance of the NAEX model such as 

y(k) =^ 0 .3a (k-D-^O .ly (k-l)-O . 02u (k-1) u (k-l)-O ,0^Q (k'2) u (k-1) 
30 '0 . QGy (k-l) u (k-3) ~0 .OSy (k-2} y (k"3 ) 

may be obtained from the general form (1) and (2) with 

12 
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coi(l)=Q.3, cio(l)=0.7, co2(l/I)==-0.02, cq2 (2, 1) ^-0 . 04, 
(1. 3)="0.06> C20 (2.3)^-0.08, else Cpg^( . ) =0 

5 Simplified desigrxs will be considered below where the 

NAKX model with only input nonlinesrlties la employed. 
However/ it will be a.ppraeia,ted by one skilled in the art 
that the present invention is not limited to use in relation 
to only input nonlinearities . The present invention can 

10 equally well be used in circumstances of both uon-linea.r 
outputs and non-linear inputs and outputs . Equally^ the 
invention is not restricted to realisation as a NARX model. 
The invention may be realised using many alternative model 
forms either in discrete time or continuous time* Models 

15 such as the Hammerstein and Wiener model, or continuous time 
models / for example, a nonlinear differential equation model 
could be used or any other model, including discrete or 
continuous spatial models / that can be mapped into the 
frequency domain. However, for each of the models the main 

20 design principle is the same. 

The NARX model with only input nonlinear! ties is given 
by equation (1) where 



2c,,(l,)y{k™X,) + 2^oi{^i)^<^-V 
1.1=1 i^-t 



for n^2 
for n=l 



(3) 



25 



30 



The structure ot the NHRX model CI) and (3) is defined 

by the values of Kjj, n^l ^ , . . , and, for each n (an 

integer between 1 and N inclusive) , involves tiarms of the 
form 



li — l/***/Kj^/ i — l/«,»>n/ 



13 



SUBSTITUTE SHEET (RULE 26) 



30-08-00 11:54 FronrHARRISON GODDARD FOOTE 



+441U2730312 



T-663 P. 21/99 Job-687 



WO 99/45644 PCT/GB99/00550 

when n > 2 and 

5 when n— 1 in the model. 

The parameter N in the uon-linear inodel ( 1 ) and { 3 ) is 
associated with the realisability of the model required to 
give effect to the energy traua formation. The ability to be 
10 able to realise the energy transformation is determined from 
the relationship between the input and output: frequencies or 
frequency ranges of non-linear systems. 

The structure parameters K^, n = l/--..N, are associated 
15 with the extent to which specific design requirements such 
as the magnitude and/or phase of the output spectrum over 
the required output frequency band can be satisfied. These 
parameters are iteratively determined as part of the design. 
The model could initially be assumed to be of a simple form 

20 in terms of these paxaxaetera . Howexrear, if the initial 

choice of parameters does not produce a satisfactory design 
the parameters are progressively or gradually revised 
according to the energy transfer effect of the resulting 
non-linear system* 

25 

For systems described by the HABX model (1) and (3), 
ttie relationship between the input and output frequencies or 
frequency ranges is given by 



30 f,= f.,UV. 

Where f^ denotes the range of frequencies of the output, and 
fy^ and ^ denote the ranges of frequencies produced by the 
Nth-order and ( ) th-ordez: noniineaxrit lea f and 

35 



14 



SUBSTITUTE SHEET (RULE 26) 



30-08-00 11:55 From: HARRISON GODDARD FOOTE 



+441U2730312 



T-663 P. 22/99 Job 



WO 99/45644 



pCTyGB99/00550 



when 
when 



lic=0 



nb 


na 




<1 


(a + b) 


(a + b) 




lib 


aa 






(a + b) 


(a + b) 






n = N 


and 



N-1 



where [.] relates to or means take the integer part, 

-J^ 1 + 1 
(a + b) J 

=[na-k(a + b) ,nb-k(a + b)l for k = 0 i' - 1, 
I^. =[0,nb-i'(a + b),], 
and the frequencies of the signal to be processed are in the 
range defined by the interval [a,b] . 

Given [a^Jb] and the required output frequency range fy/ 
the smallest N for the NARX model (1) and (3) which can 



domain 



bring about the spec=i£ifed frequency 
transformation can be determined from equation 4- 



energy 



13 



20 



25 



(iv) Map the NARX model with the structure given in (iii) 
into the frequency domain to yield the frequency domain 
description. The frequency domain description is given in 
terms of the Generalised Frequency Response Functions 
(GFRFs), H„(jw,,---,jw„) , n = l,...,N, which, after this mapping, 
are specified in terms of time or spatial domain model 
parameters . 

The mapping of the NARX model (1) and (3) between the 
time or spatial domain and the frequency domain is given by 



H.(jWi,--vjw„)= 



n = 1/...,N 



(5) 



15 
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Therefore the frequency domain proper'cies of the system can 
be completely defined in tsrins of the parameters ( . ) of 

the time or spatial domain description of the system. 

5 (v) The output frequency response o^ tXie non- linear system 
(1) and (3) is given by 

n-1 



10 



where 



with 



J( -) do. 



15 denoting an integration over the nth-diiuensional hyper-plane 

Based on this relationship, the parameters in 
(jwi^ . . . ^ jWf^} . . . ^N^ which^ due to the mapping 

20 performed in (iv) , are the same parameters as those in the 
time or spat:ial domain model are determined. This step 
enables the shape of the output frequency spectrum Y(j^) to 
be defined which in turn ensures that the Spectrum 
approaches the specified output frequency spectrum as 
25 closely as possible. 

Different design specifications can lead to different 
implementations of corresponding designs. 



(v.l) Firstly, given knowledge of the input spectrum 

30 U(jw) and the required output spectrum (jw) . Substituting 
(5) and (7) into (6) yields 
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J exp[-j (w,lL+/-'+w„ljinu(jwJdcr„ 



Equation (8) enables the parameters associated with the 
time or spatial domain NARX models 

and 

^10 (^\^ ' -^1=^/ • • • ' ^1, 

to be derenuincsd as follows to implemsnt the required 
design : 

1) Based upon the equation 
^ l/Vn ^ 



15 



determine par ante ters 

using least sq|;ua.ras routine to make the right hand side of 
equation (9) approach the specified output spectrum as 
20 closely possible • 

TThe first term on the right h^nd side of equation (8) 
1 



Scio(l^) exp (-jV7l^)| 

25 

is omitted from (9) . The omitted te2nn represents linear 
output t^rms in the time or spatial domain realisation of 

17 



SUBSTITUTE SHEET (RULE 26) 



30-08-00 11:55 Froni:HARRI SON GODDARD FOOTE 



+441142730312 



T-663 p. 25/99 Job-S87 



WO 99/45644 PCT/GB99/00550 

the nonlinear system and these may not be needed to achieve 
the design at this step, hence they are omitted from (9) . 

2) In order ro augment the performance of a filter designed 
5 as above, it may be desirable to also design, a suitable 
linear filter H(j^) to improve the approximation to Y^ijw) 
obtained in 1) above such that 

can achieve a better approximation to Y* {jw) , As part of 
10 this linear design the parameters c:xo(-^l^^ -^1 ^ly-^./i^i^ 
which ware omirred from (8) to get (9)/ can be obtained from 
the parameters of the linear filter. 

Design 1 hereafter illustrates the design of a non- 
15 linear system using this first case. 

(V,2) Secondly, given the input spectrum and a, 

specified bound for the magnitude of the required output 
spectrum; Y^"*" (p^) , 

20 A bound ( w) for the magnitude of Che output spectrum 

Y{jw) at the NARX model (1) And (3) can be expressed as 

Y^w> = 2;;— ^|HJjw,,--vjwj|°lu|*-"*|n(jw)| (10) 

25 according to the result in Billings, S.A, and Lang, Zi- 
Qiang, 1996/ A bound for the magnitude characteristics of 
nonlinear output frequency response functions, Part 1: 

18 
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Analysis ana camputation . Int. J. Cont:roi, Vol. 65^ pp309- 
328, where 

|Dl*'-*lU(jw)| 

5 denotes the n-dimensional convolution integration for the 
roagnitude charactezristic of the input specrruia and 

represents a bound for the GFRFs magnitude 

10 with W]^ / - , w Pv^n satisfying the constraint wi+, - - - , -^Wj^ = 
For the HARX model (1) and (3), 

can be e^raluated. as follows 



15 |hJ jw,,-", jwj|^ = 



Jin 



1" Xci(jCl^)exp ("jwl,) 
Combining (10) and (11) yields 
1 



;i2) 



20 



where 

25 Equation (11) enables 

and cio ^-^li' Ii-l^-.-^iCl to be determined for shaping the 
bound {w) for Y(jw-) in order to make this bound approach 
YS*(t^), The procedure which e^a.n be used to achieve this is 
30 as below. 
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1) Based upon the equation 

= CJUH— *iU(jw)j (13) 

use a Isast squares routine to make the right hand side of 
S th^ above equation approach Y^*(v) especially over the 
frequencies or frequency range to which the required energy 
transformation needs to be iinplemented. The coefficients 
n-lr-'r^ f in equation (13) must be constrained to be 
positive since is^ the result of the sujcimation of the 

io modulus of the coefficients <^0n (^1^ ^ ^ ^ r 

The first tei-m cn the right har^d sid^ of equation (12) 
_1 

l-S<=io<^i> {-jwlj 

15 is omitted from (13) . The omitted term represents linear 
output terms in the time or spatial domain realisation. 
These omitted terms may not be needed to achieve the design 
at this step, hence they are omitted from (13) , 

20 

2) If necessary, the approximation to Y^*(p7) in 1) above 



can be supplemented using a linear filter with a maguitude 
25 characteristic lH(jw)l such that 

^ ^/t; ^ • V 

provides a better approximation to the specified bound and 
as a result cxo<^l)/ Ii=l, * . - , K^, which were omitted from 
30 (12) to get (13) , can be obtained from the linear filter 
parameters . 

20 
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Design 2 hereafter illustrates the detailed procedure 
of the above and several examples of the design. 



(V.3) Thirdly, there are many practical situations which 

5 should be dealt with on an individual basis. There follows 
examples of two such situations . 

(a) Referring to figure 3 there is shown an input 
signal having a spectrtim U(jw) 300 which comprises .a portion 
10 302 of interest, between frequencies a and e, and an portion 
304, between frequencies e and b, ro be translated to 
another frequency range 306 defined by frequencies f and 2b 
whilft retaining the portion of interest 302 within the 
antput signal 30S . 

15 

In order to realise the energy transformation shown in 
figure 3, a non-^linear system 400 can be constructed as 
illustrated in figure 4 which comprises means 402^ 404 and 
406 for implementing the following components 

20 

H^(jw) 402 and H(jw) can b© iiapleiuented readily using 

classical linear band pass filters, while Hzij^lf j>^2^ 404 
can be constructed in the time or spatial domain as shown in 
25 the equation below 

with the parameters, c^jdj^^lg) , li = 1/---/K,,l3 = being 
determined to produce the signal y2(^> with a required 
30 frequency characteristic. Each of the components 402, 404 

21 
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and 406 have corresponding frequency responses 408, 410 and 
412. Consequently^ the whole non-linear system can be 
realised as a non-linear time or spatial domain filter as 

y(ic) = Hi{q"')u(k) + H{q-') £ c^^C Vlj)];!^ (k - 
where 

are the backward shift operator descriptions of the linear 
filters vhich have the frequency response functions Hj_ (jw) 
and H(jw)f respectively. 

15 The equarion foj: y(k) can be fur-cher wxrinT^en 

Hip(q"')H^(q'')y(k) = H^(q^')H^(q-')u(k>+H„(q"')H^{q"h £ Cp^d.^uH^ (k- IJ 

which is clearly a NARX model that can be described by 
equations (1) and (3) . 

20 

Although the above shows energy transformation from a 
frequency band to a higher frequency band, energy can 
equally well be transferred from one frequency band to a 
lower frequency band. 

25 (b) If the objective of the energy transformation and 

hence the desired non-linear $ystem is only to distribute 
energy of the signal ro be processed over a wider frequency 
band without energy amplification then a simpler model may 
be sufficient. Referring to figure 5, there is shown 

30 schematically an input signal SOO in the frequency domain 
comprising energy between frequencies a and b. The desired 
output frequency spectrum 5 02 comprises two portions, a 
lower frequency portion 504 and a higher frequency portion 

22 
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506, It: will be appreciated that the energy of th-e input 
signal SOO is to be spread ov^r the two frequency ranges 50 4 
and 506. 



can be designed to redistribute energy of the original 
signal u (k) , with fr^Quency components over the frequency 
10 band [a,b], to the new ranges [Q^b-a] and [2a^2b] without 
energy aitiplif icatiou provided an appropriate a is selected 
in the design process described above. 

(vi> The final step in the design process is to 

15 materially realise^ that is, to physically realise the 
designed filter using appropriate software or hardware or 
coirtbination of software and hardware. 

Although the present invention has been described above 
20 with reference to a NARX model, it will be appreciated that 
the present invention is not liiuiiied thereto. The method 
and realisation of non-linea.r systems having pre- 
dete2rminabl^ frequency or energy transfer characteristics 
can equally well be utilised using other forms of 
25 descriptions of non-linear systems or models where there 
exists a requirement to transform or transfer energy at one 
frequency or band of frequencies to energy at another 
frequency or band of frequencies. 

30 The present invention can be used to realise energy 

transformations over frequencies using a nonlinear system or 



5 



A quadratic filter^ for example. 
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to modify energy rraixsf ormarion of &n existing linea.jr or 
non-linear -system. 



Furthermore, the present invention can be applied or 
5 utilised In Che field of design and r^ali:5etion of 
electronic eirouits or filters. The energy in tine input 
signal at a particular frequency or band of frequencies is 
transferred to desired frequency bands. Similarly, in 
mechanical systems^ t.)ie addition of non^linear mechanisms 
ID could transfer the energy of a vit>ration at an undesirable 
frequency to some other frequency. The present invention 
may also find application in the field of fluid mechanics, 
for example/ in the effects of flaw around objects (e,g, oil 
platform legs), noise in ducting and pipe flow systems. 

15 

Alternatively, the modifications which are required to 
be effected to a known linear or non-linear system, for 
example, a mechanical system, in order to bring about a 
particular frequency distribution of energy can be 
20 determined using the present invention and then the linear 
or non-linear system can be so modified. 

25 XX -1 DBSJGW 1 

Design No,l shows an example of how the energy of an 
input signal having predetermined frequency components can 
be transferred to other frequencies. 

30 

Consider digitally filtering an input signal 
u(t)=:t cost 4-cos2t, (14) 

24 
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wherein the sampling period is T^^ 1/100 s- 

The first step is to determine the frequency spectrum 
o£ the signal to be processed. The frequency spectrum of 
5 the signal to be processed comprises input frequencies 
w^i - 1 and =^2 

The second step is to specify the output frequency 
characteristics of the filter. Two different filtering 
10 problems will be considered ror this example. 

The specification for the first problem is to transfer 
the energy in u(t) to Tihe outpuc frequency 

specification for the second problem is to trxransfer the 
15 energy in u(t) to the output frequency =4 . 

Since the sampling period is T = 1/100 s, rhe digital 
input to the filter is 
1 2 

u(k) = cosY^k + cosY^k for k = 0,l,-.* , (15) 

20 the normalised input frequencies are therefore 

w^, =1/100 and w^^ =2/100, 

and the required normalised output frequencies are 

w^^=0 for the first filtering problem and w^^ ==4/100 for the 

second filtering problem. 



25 



For this example, the output frequencies produced by 
the nth-order nonlinearity of a non-linear filter are 
distributed uniformly in 
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10 



IS 



20 



25 



£ _ 



When 



+ b) 



( a + b) 

L a 



< 1 



(16) 



> 1 



Ml, when 

X ] means ro take the integer part 



1 2 

where a = , b=" , [ 

100 100 

of X, and 



1 = 



na 



(a + b) 



I. =[na- j (a + b) , nb-j (a + b)] for j = 0/ ,*,/i*-l , 
I.. -[o,nb-r(a + b)], 

with the difference between any two neighbouring frequencies 

2 11 

being , 

^ 100 100 100 

For n=l^ it can be evaluated from the above equations 

(17) 



that 



lOO'lOOj 



and the corresponding output frequencies are 

_1 2_ 

100 ' 100 

For n~2, it can be similarly obtained that 



2 


4 










_100' 


100_ 





0, 



100 100 



lOO'lOO 



0, 



100 



and the corresponding output frequencies are 
12 3 4 



0. 



lOO'lOO'lOO'lOO 



Therefore, a NARX model of u.p to second order 
noniinearity will be sufficient to realise the energy 
transformations required by the filtering example thereby 
addressing the third step in the design process. 

According to t:he above analysis, select a ^^ARX inodel of 
the form 



(18) 
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as the basic structure for the non-linear filter and take 
K^r^l for simplicity. The parameters K^, c^o(l) / — /Clq(KJ and 
c^^ik^r^:,) f k^^O.l and k^^Q^l, then need to be detenaiued to 
complerely specify the design. 



10 



1 ,i 



15 



20 



In or-der to derive th^ procedures for determining these 
parameters, consider the frequency domain characteristics of 
the filter model (18), According to J, C.Peyton Jones and 
S.A.Billings, Recursive Algorithm for Computing the 
Frequency Response of a Class of Non-^linear Difference 
Equarion Models, Int. J. Cont:rol, 1989, vol.50, No, 5, 1925- 
1940, the generalised frequency response functions of this 
filter model ^r^ 



H^( jw>-0 



22^02^ Kr^nJ exp[- j(w,k, + Wak^)] 



1 - i Ci^^( k,) exp[- j(w , + W2)k, ] 



(L9) 



The output frequency response of ths filter (18) 

under the input 

1 2 ^ A(w.) 

u(k) = cosY^^ + cos - ^ 



i — k= y 

100 ..iri 



(20) 



with K=2, A{w.) = l, w^=^QQ, 



:±1,±2> can be 



25 



written as 



?(jw) = SYn<3w)-Y2(jw)- 



^ -J 



30 



for W'> 0 

for w = 0 
(21) 
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wYiere 2,-1.-^1 ,+2) . Y{jw) is related to the output 

spectrum Y(jw) by 

r2Y{ jw) for w >0 
^^^^^^^lYCjw) for w-0 

and. similarly, Y,(jw) is related to the nth-order output 
5 spectrum Y^(jw) by 

r2Y^(jw) for w>0 
j^) ^ 1 Y ( jw) for w = 0 
This is a" specific form of the general expression for output 
friquencj responses of non^linear systems which are g.vcr. by 

^ . . ' (22) 

Y(jw) = 2^Yjjw) 

10 where 
with 

de^^o^lng an integration ovar the nth-dimensional hyper-pLane 
15 v/,+,-,+w, = w and N being the maximum order of the dominant 
system nonlinearities . 
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PCT/GB99/00550 



- — i;i;c„(k„k3)^ za^^,^) a(w^) expf- j (w, k, k^)] 

l-ic,„(k,)exp[-jwk,] — 



I I 



= H(3W)l;Ec,,(k,,k,)f (w,k„k2) for w>0 

and 

Y(jw) = H(jw)tt=o.(lc.<k,)^^^^^f^ for w^^O 
where 



H( jw)= 



l-S^io^J^L^ exp[-jw ki] 



(24) 



(25) 



(26) 



10 



f (w,k,,k,)=-^ SA(Wi,)A(w,^)exp[-j (Wik,+Wikj)] 



Moreover denoring 
i i 

S Z^oz( ^i'J^2) f (w^k.^k^) for w > Q 



Y ( jw)=^ 



1* I 



, £(w.k,,k,) 



for w=0 



gives 

15 y (jw)-H(jw) Y(jw) . w^O 



(27) 



(28) 



(29) 



In view of ChS fact that HCjw) is tUe fcecjuency 
response function of a classical linear filter and Y{jw) 
a linear function oC the filter pararo*ters Cq^{}^^,)c^) , 
20 k, = 0,1 and k^^O,!, the procedures for determining the 
parameters of the non-linear filter with the given design 
requirements and the structure (18) are given as : 
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(a) From the design requirements, determine the desired 
output frequency characteristic Y^(jw) and choose the 
parameters CQ^Ck^.^) , ki^=Q,l and k2=0,l, appropriately to 
make Y(jw) a.ppro>:imate Y*(jw) as well as possible, 

5 

(fa) Examine the filtering effect of 

yW^Y,T.^02(KK)^i^'K)^i^-K) f (30) 

where CQjfkj^k^) / k^^^O,,! , k^^O^l, are the results obtained 
in the procedure (a) ♦ If the effect is accepcable, then 
10 choose H(jw)=l so tliar the filter parameters 

otherwise design the classical linear filter K(jw) to make 
H(jw)Y(jw) satisfy the requirements for the output frequency 
characteristics and obtain the filter parameters 
15 K^/ c^^(l) ;C^q(2) /,..,Cio(Ki) at the same time. 



Therefore, in ordar to address the firs^; problem to 
transfer energ-y from u(t), w^j_ = 1 and w^^ - 2 , to the output 
frequency w^^, = 0 , take 

20 r(jO) = l, Y*(jwJ-0, w.=Y^, i-1,2,3,4. 

and for the second problem to transfer energy from u(t:), 
w^3_=l and w^2=2, to the output frequency w^^=4, take 

Y'(j^)-1, Y'(jwJ = 0, w,=~, 1=1,2,3,0- 

25 The filter parameters o^^^^xf'^z) ' k^=0,l , k^^Cl, can Chen 
be determined through the group equations 
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Re Y (jO)]=2.ZCo2(^i'1^2^ 0 



Re 



1 1 



10 0 J v:i=<lfe3-0 



Im 



"0161=0 



(31) 



fay the. least square's method. In (31) •) and f'( .) 

represent the real and imaginary parts of f{.) ■ 



Rewrite (31) as 



(32) 
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where 
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Re[Y*(jO)],---,Re 



j-- — ) 
^100 



,lxa[r(jO)]r-,lm\Y'{j^) 



f ^( Q , 0 , 0 ) 



f (0. 0, 1) 



.^(^.0.0). ^^W^'^^ 



f ^( 0 , 0 , 0 1 



f V 0 , 0 , 1 i 



^''iT~^^^^^' f^TXTT'l'l)^ f'^TTT'O.l), 



100 

2 



IQQ 



100 



2 c o s 



100 
2 

10 0 
3 

10 0 
4 

10 0 



2 COS — TTT 

— 

''2 1 

COS _ _ _ + G O S 



0 *- 2 3 1 n 



— Sin 



0 - 2 s i n 



X 0 0 
2 

10 0 
3 



100 ~~ 100 
1 

^ ^ s 

10 0 

I 2 
cos _ , _ + c o s . ■ ■- 
10 0 10 0 

2 

COS 

10 0 
i 

2 ^ ^ Tirr 
1 0 Q 

2 

1 2 

= ^ ^ „ „ ^ "Si:: 

1 0 Q 100 
L 

— Sin — ■ 

LOO 



10 0 
1 C 0 



sin 



10 0 



X 0 0 



(33) 



(34) 



(35) 



X 0 0 

-sin 

Notice that 9 can be written in the form of (34) due to the 
fact that 

10 f {w,0,l)=^ f (w,l,0) , (36> 

Therefore, the filter parameters 9 for the first 
filtering problem to transfer energy from u(t)y 

~ ^ ^az = 2 / "the output frequency w^^ - 0 can be 

X5 Qbtalned as 



4 ^(x^r'x^Yi 

where 



(37) 



1,0, "V 0,0, --v 0 



(38) 
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and. the parameters d for the second filtering problem ro 
transfer energy from u(t), w^i=l and w^j=2, to the output 
frequency =4 can be obtained as 



4 =(X'X)"'X-Y2 
5 where 

4 s 

Y, = 0,"-, 0 1 O,'--, 0 



(39) 
(40) 







1137 




C^zd, 1) 




1187 


(41) 


0^2(0, 1)+C^2(1,0)_ 




-2373.8^ 








"-1206.2"^ 




C-02{1, 1) 




-1206.2 


(42) 


_C'02(0/1) + C'02(1/0)_ 




. 2413.2 ^ 





The results of computing (37) and (39) are 
A. 

and 

10 4 = 

respectively. 

Under the input of (20), the power spectral densities 
of the input and output of the non- linear filter 

15 yi(k) = 1; 2c^o2(k:,,k2) u (k - Jc^) u (k-k^) 

= cS2(0,0)uHk) + £'o.(ia)u"(lc - 1) + [£'32(1,0) + c^2(0,l)]u(k - l)u(k) 

(43) 

which is initially designed to address the first filtering 
problem above, are shown in figure 6, and the power spectral 
20 densities of the input and output of the non-linear filter 
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i i 

= c%^(0, 0)u^(k)H-c'o2(l/ l)u'(k-l) + [c%2(l,0) + c'o2{0,l)]u(k"l)u(k) 

(44) 

which is iuitiaiiy designed to address the second filtering 
5 problem above^ are shown in figure 7 . 

Further improvemem: to the performances of the filters 
(43) and (44) may be possible. Therefore^ a linear filter, 
H(jw)^ is designed in order ro further imprcve the filter 
10 perf o2:mances - 

To improve the performance of the filter (43) for the 
first filtering problem, H(jw) is designed to be a fifth 
order low-pass type 1 Chebyshev filter with cut off 
15 frequency 0*5 rad/sec and 0.5 dB of ripple in the pass-band 
to ensure a satisfactory frequency response over the pa^s 
band. The result is 



b'(l) + b'(2) z"'+'"-+b'(6) 2 



-5 



H, ( jw) = : -T T 

- -i^Q-L^ Q-^^'^^ + Q-QQ^^^""^ 4-Q.13I22~^+Q.1652g"^ ^0,Q933z^' + 0.01592 
l + 4,994l2~^-9*S765z'^-h9. 9648z"^-4. 9766z^' +0 .99422"^ 

20 ( 45) 

To improve the parformaneet of the filter (44) for the 

second problem^ HCjw) is designed to be a fifth order high^ 

pass type 1 Chebyshev filter with cut off frequency 3*9 

rad/sec and 0.5 dB of ripple in the pass band for the same 

25 purpose as in the first problem case. The result is 

V ( 1) + b^ 2 ) -h- - -+b^{ 6 ) 2"' 
' l"a2{2)z"' a'[6)z-* 



0 , 9218-4. 60882"' + 9 . 21752"^ --9 . 2175z~^ -h 4 . 60 8 8 a"** - 0 - 9218 2"^ 



l + 4.3381z"^^9.3635z*^+9.0613s'^-4 .38 4 6z"^-i-0.8 43 6z~^ 



(4S) 



30 The purpose of the additional linear filter is to 

attenuate unwanted frequency components in the outputs of 
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the filters (43) and (44) for the above two filtering 
problems respectively zo make the output of the additional 
filter satisfy the corresponding design requirement. 

5 It will be appreciated that the filter parameters 

K^, o^o^l) ,c^^[K^) f associaT:ed with the expression 

|l- ^CifldCi) exp[- jwkjj in H(jw) can be obtained as a result 

of dividing the denorainator of H{jw) by the numerator of 
H(jw)* The specific H(jw) is given by K^(jw) and H^Cjw) in 
10 equations (45) and (46) for the two filtering probleicis 
respectively . 



The general description for the nonlinear filters 
15 designed as above is 



20 



ty(k)-H(q"')y (k) 



y(l<)-£c,o(ic,)y (k-k,)- l-2l<^,ciK 
whore 

HCq"^} = l-2c,,{k,)q-*M 



(47) 



(43) 



and q" , and q^ denote the backward shift operators . 

Another expression in the time domain for (47) is 



25 Therefore, eiribodiments of the non-linear filters which have 
been designed can be realised in a manner as shown in figure 
B or more specifically as shown figure 9. 



It urill be apprecxattfid that In figure 8 ttie two 
30 components 800 and 8 02 are represented by 
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and 



1-a (2) q^'— a (6)q"' 



For the first filtering problem in this example 
With Co2{ki,k2) 's given fay (41) and 



. bNl) + b^(2)q"^-i""+b^(6)q- 
'''^ ' l-a^{2)q'^-"-a^(6>q-^ 



with b^<i) 's and a^i) 's giveii by (451 - 

J 

For the second filtering problem in this example 

with c^2(k,,kj 's given by (42) and 



( 1) + b" ( 2 ) q-^ +- - --f b' ( 6 ) q"^ 



l-a^2)q-^-— a=(6)q- 
with b^(i) 's and a^(i> 's given by (4^) . 

Referring to figures 10 and 11, there is shown the 
filtering effects of the filters designed to address the two 
filtering problems. Figure 10 shows the power spectrum 
densities of the output and inpur of the nonlinear filter 
which is finally obtained for the first filtering problem. 
20 Figure 11 shows the power spectrum densities of the output 
^nd input of the nonlinear filter which is finally obtained 
for the second filtering problem. The effects of the 
filtering in the time domain of the two filters are shown in 
figures 12 and 13. All the responses indicate that the 
25 filters sub sta^nti ally satisfy the design requirements. 
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XI .2 DESIGN 2 

Design No. 2 illustrates the detailed procedure and 
several examples for designing a non-linear system for 
transferring the energy o£ a. signal from a firar 
predeterminablc frequency or rang© of frequencies to a 
second predeterminable frequency or range of frequencies 
such thac the output frequency response of the nonlinear 
system so designed is within specified bounds. 



II . 2 . 1 Oetailed Proeeduare 

;= (1) Given u(c), the signal to be processed in the time or 

lJ spatial domain, the frequency band [c,d] to which the 

i IS energy of u(t) is to be transferred, and the user 

^ specified bound Y®'(w) for the output spectrum Y(jw) over 

[G,d] for the design. 

^ (2) Sample the time or spatial domain signal u(t) with 

3 20 sampling interval T to yield a discrete series (u(k)}and 

^ perform a Fast Fourier Transform (FFT) on the series to 

compute the spectrum U(jw) of u(r) as 

<^^)'-^m ^-ii-) 

where U^[j (■)] is the result of the FFT operating on {u(k)} 
25 and M is the length of the data used to perform the FFT. 

M is taken as an even number for convenience, 
(3) Evaluate the range [a,b] of frequencies in u(t) as 



2fe __ 27T 

MT ' ^ ~ MT 
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where is an integer such that 



PCT/GB99/00550 



v-'mt 





\ m J 


an integer such 


>0.05, 





<0.05 for i e^dt, +1) r-'/ — r 

2j 



(4) The relation.ship between the bound of -che output 
spectrum Y*(w), the coefficients of the NARX model 



2lci3(li)y<k-lJ+£co^(l,)u(k-i,) for n=:l 



and the spectrum U(jw) is given by 



V = 



where 



i-i, 1 



are-^arameters associated with the NARX model parameters 
ju|*-"*jij ( jw)l denotes the n-dimensional convolution 

a 

integration for the magnitude |u(jw}jof the spectrum 
U{^w)f which is defined by 
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and No=l when the NARX model involves nonlinear terms 
from order 1 to N . 



!0 



15 



20 



25 



Based upon this expression, the structure 
parameters N and and the NARX model parameters are 
determined as below. 



(i) Evaluate 



[Jl^ when 

1=1 



wnen 



(a + b) 
nb 



Li-o 



na 



(a^b) 
na 



<1 



> 1 



where [ x ] denotes the xntsger part of k, 

, . r 

L (a + b) ^ 

I, ==[na-i{a + b) ,nb-i{a+b)], for i = 0,-vi'-l 

I^. =[o,nb'i'(a-f b)] 
for n"l,2,... until a value of n is reached such that 
part of the specified output frequency range [c,d] 
fails into f^ . This value cf n is used as the value 
of N, . 

(ii) Evaluate 

^or n=^2,3,... until a value of n is reached such 
that the frequency rang-© [^fd] falls completely 
within the correspcinding fj. . This v-alu.e of n is 
taken as the value of N . 

(iii) Calculate 

M*'-^*|u(jw)| 
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to yisld 

lul^*-'^ju(j2;r i/MT)!, for i = 0,--vM/2 



10 



15 



20 



using the algorithm 



2is 



0(i) = 



for n = No,**->N, where Conv(.) denotes 



convolution operation and U(,) and U ( .) rftpreaent 
the intermediate results or this algorithm* 

(iv) Based on the equations 



where i, = round 



CMT 



Z3Z 



TdHTl 

; id = round ~z — * ^rid round { x ) 

means to take the integ-er nearest to use a least 
squares routine to compute 

under the constraint that the results must -be 
positive, and then select the MARX jnodel paraiaeters 

1: =l,-^K„*-.l„-l."vK, 
for n=^No,-vN 

•^der the constraints on the summation of the 
modulus of the coefficients given by 

(v) If necessary, design a classical linear filter/ for 
example/ a band pass filter which ideally allows the 
frequency respons'e to be unity over the frequency 
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band [c,d] and zero beyond to yisld the linear 
frequency characteristic 

1 

J 

and therefore . determine the parameters associated 
3 with K^, c^qII,}, l^ = l,-",Ki, Otherwise, all of the 

parameters of c^^l.) can be taken as zero to yield a 
model having np regression terms associated with the 
output , 

(5) Construct a NARX model as shown in figure (14) using the 
10 results obtained in the above (iv) and (v) . The 

nonlinear system illustrated in figure 14 comprises a 
nonlinsar part 1400 and a linear part 14Q2. It will be 
appreciated that m figure 14 

i ^c^(i:^"vla)flu(lc-l.)=i;-E <^{l,r-X)th(^-k) 

15 

11,2.2 Three Specific Examples 
Example 1 

20 This example illustrates a further implementation or 

design of a nonlinear system following the above detailed 
procedure . 
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(1) The signal to be processed is given by 

1 sin at -sin B t 
,,,,,2M.- 

wirh a = 3, 3, /?=1, andM^,=1.6, The frequency spectrum 
U(jw) of the signal is shown in figure 15 indicating 

5 thai: the real input frequency range is [l,3.3]. The 

requirement for the design is to transfer energy of the 
original signal to the frequency band [c/d]= [5 . 6, 7 . 6] 
with the bound ^ on the output spectrum magnitude 
specified to be Y°*(w) = l*6 over this frequency band. 

10 \2). Sample uit) with sampling interval T = 0.01sec to produce 

1 sin ^xkT-sin ;^ kT 1 sin 3 .3i^ 0, CXk-sin 0,01 k 

u{k) = 2M^ = 2xl.6x— X — — 

^27r kT 2jr 0,01k 

k = -1993,'*v0,'-v2000 
and perform a Fast Fourier Transform (FFT) on this 
15 series (M=4000) to compute 

u/j — ll-u/j-^ll l--1999r".0r-*,2000 
and -hen to yield 

r/j— il=[/j .^^^ ^ 'lUo^oiu/j-^ij 1=^1999,", 

J V4Q00x0,01 J V4000 J 

the result of whxch, in the nonnegative frequency range, 
20 is shown in figure 16. 

^otxce the difference between the real spectrum of 
u(t) in figure 15 and the computed spectrum in figure 
16. The differences are due to the errors caused by the 
25 FFT operation. The design should and will be performed 

based upon the computed spectrum to lead to more 
practical results. 
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(3) Evaluate the frequency range [d,b] of u{t) from the 
computed spectrum giving 

a = 0.6283 b = 3.7699 

(4) Design the system structure and parameters 

5 {i) Determination of , 

Cleariy> the output frequency range contributed by 
the linear part when the input frequencies are within 

[a,b] = [0.62a3, 3.7699] is 

=[a,b] = [0.6283, 3.7699] 
10 The frequency range f^, produced by the second 

order-" nonlineariry in this case is obtained as 
follows , 

since n-2, 



nb 

15 (a+b) 



20 



na 

(a+b) 

and 



2x3.7699 
'(3.7699+0.6283) 



2x0,6223 



(3,7699+0.6282) . 
1.7143-1 0.2857 ]=1,7143-0>1 



1 - 



na 



2x0.6283 



_(a+b)J I 3.7699+0.6283 J 



^1==[ 0,2857 ] + I-0+l^l 



So 



4se,ii|j0^{a+bm2x0.62S3, 2x3.769^0, 2x3. 7699-(0. 6283+3.7699)] 
=:[l.2566. 7.S399jJ[0, 3.1416] = [0, 7.5393] 
fy^ thereby obtained contaiin:s pare of the specified 
output frequency range [c/d]-[5.6, 7.6]. So, is 



determined to be Nj = 2 . 
(ii) Determination of N , 
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PCT/GB<^9/00550 
Evaluating f, =t^[jf^_^^^ tor ri-2 yields 

^Uz=4U^.-[Q/ 7.5398jU[0,S283, 3.7699] = [o, 7.5398] 

To evaluate t,\^, ^ , calculate f,^ first. In 

this case 

^ ~ 1 



_nb 



1 = 



na 
(a+b) 



3x3.7S99 



3x0,5283 



na 



{a+b) 



(3.7659+0.6283)1 (3.7699+0,6283) 

=2,5714-[ 0.4286 ]-=2,5714-0>l 

3x0.6283 
. 3.7699+0,6283 



+ 1 = [ 0.4^86 ] + l = 0+i = l 



So 



="U^. =IoU^i =[na,nb]U[0,nb-{a+b)] 

= [3x0,6283, 3x3.7699U[0/ 3x3, 7699-(0. 6233 + 3. 7699)] 
=[l,884S, 11.3097]tj[0. 6.9115]=[0, 11.3097] 
Therefore 

^U=4,Ut.-[0. 7.5398]U[0. U.3097] = (0, 11.3097] 

^vln^i thereby obtained includes the whole spcscified 

output frequency range [5.6, 7.6], N is therefore 
determined to be N = 3, 

(iii) Calculate |oK---|a ( jw)| for n N, = 2 and n = N - 3 ; 
respectively to yield 

|t^-^U(j27 VMrj={l|^>^U(j2g j/4QOQxO.Qlj i^O,- >400(?2, n-2 ard 3 

The results are shown in figure 17. 

(iv) -Based on the (i^-i^+l) equations 

Y^*(/r i/2000x0.01) = 1.6=:-^C2 |ui*|u(j ;r i/2000 x 0 , 01 
+ -^^C3 |aHui*|D{j;r i/2000x0.0l)| 



with 
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= i:ound[cMT/2;r] ^ E:ound[5 .6x4000x0.0 1/2 ?<:7r] = 36 
1^ = round[dMT/2;r] = rouad["7 .5x4000x0.0 1/2 xtt] = 48 
use a least squares routine under the constraint of 
nonnegarive solutions to compute O and C3 , that xs, 
5 to determine c. and C3 , under the constraints of 

CjSO and C^zO, ro nLinimises the following 
expression 

jfl.e-ic, |t:i^i:(j;fiy2000xO.Ol]l--^C3 !l|*it^%;ri/2000x0.0l| 

The results obtained are 
10 ' €2=0 and C3 =3,8367 

(v) Desxgn, oprionally, a linear Butterworth. band pass 
filter to attenuate the frequency components beyond 
the frequency band [c, d] = [5.6, 7.6] to yield the 
linear frequency characteristic 



15 



10-5fo,098S-0.1972q"- 4-0.0986q"') . 

!^ . ' TH-^ 

l-3.9633q"^ +5 . 8988a'^ -3 . 9076q-^ +0 . 9721q" * 
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(5) Constroict a NARX m^del as shown in figure 14 with 

N, = 2. N = 3, K,=K,-1. c,,(ia) = 0, c.^l 1, 1 . D = ^ - 8367 , 

that is, 

L I 



- =c^aA)^^k-lH=b,aA,lla^(!c-l)=0xa=(k-l)+3.8367u'(k-l)=3.a367uP(k-l) 

and 

l[r^( Q.0986-0.19"72q---r0 .OSSSq") 
U-i=;oU,)q" = i-3.9633q-'+5.8986q-'-3.9Q76q'+0.9721q^ 

which determines the parameters associated with the NAHX 
10 model: parameters K, and c.odi) , li=lr-vKi. 

The output frequency response under the given input is 
shown in figure 18 indicating T:hat the energy has been 
transferred to the specified frequency band [c, d] = [S.6, 7.6] 
15 with the magnitude of the response below the specified bound 
1.6. - • 

The frequency response of the above design is examined 
for other input signals below. 
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In -th^ first: cas&, consider 



u(t)== ^ 



2co5 cosbt-Gosa t 

2 



where and are defined as above. The frequency 

spectrum of the signal is shown m figure 19. The frequency 
5 range is clearly the same as that of the signal havxng tihe 
spectrum given in figure 15 and nhe magnitude of tiie 
spectrum also satisfies zh^ condition tihat 

lu(jw)l<K, =1.S , 
This implies than the frequency response of the designed 
10 sj^stem to this u(t) should in theory also transfer energy 

into the frequency band [c, dj = [5.6, 7.6] with the output 

raagnitude frequency response being less than Y°'{w)-1,6 over 
this frequency band. Figure 20 shows thxs frequency 
response and indicates char the actual result is consistent 
15 with the theoretical predictions. 



For the second case, u(t) was taken as a random 
process with the frequency spectrum given in figure 21. The 
frequency spectrum is substantially within the frequency 

20 range [l, 3,3] with a magnitude of less ::han 1.6, Therefore, 
the same conclusion should apply for the output magnitude 
frequency response of the designed system to this random 
input. i'igure 22 shows this response and indicates rhat rhe 
energy is transferred to a new frequency band of 

25 substantially [5,6, 7»6]. Note th^it the magnitude of the 
0UT:put spectrum over this frequency band is well below the 
specified bound Y^'(w)-*1.6. This is because^ of the effect of 
the attenuation which is due to intrakernel interference of 
the nonlinear mechanism. 

30 The nonlinear filter which has been designed aBove can 

be represented by the block diagram shown in figure 23. This 
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could b© realised electronically by following the approach 
used to realise Design 1 ia figure 9. However^ for some 
applications the design may need to be realised in 
continuous tiiae- Using the bilinear transformation 
_ 1 + (t/2)s 

where T=0.01 is the sampling interval^ a"^ is th^ delay 
operator, and s is the Laplace Transform operator and 
substituting in ths discrers time design in figure 22 
provides the equivalent - continuous time system shown in 
10 figure 24. Simulating the system in figure 24 produces 
almost identical responses as the discrete "ime system 
response, of figure 18. 

The system in figure 24 could for example be realised 
15 mechanically as illustrated schematically in figure 25 where 
the cubic device is either a material which exhibits a cubic 
response or is impieniented as an actuator which takes u^(t) 
as input and produces an actuation output n^it) which is 
proportional to the cubic power of the mpui:. 

20 

One possible application of this design would be m 
vibration isolation. For example it may be required t:o 
transfer energy from an input frequency range [1^3.3] ro the 
frequency range [5.6,7*6]. The new design could be used to 
25 achieve this effect* 

Othwr much more complicated designs for vibration 
isolation can be achieved based on the present invention. 
The design procedure would be exactly as described atoove. bur 
30 the realisation would involve tine syuthesla ox dampers, 
dairying materials or actuators with the nonlinear dynamic 
characteristics specified by the designs. 

Example 2 
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Bxampl^k 2 shows an application of the design ahoT^e to 
the attenuation of signal energy over unwanted frequency 
bands using designed nonlinear effects, 

When using linear structures for the attenuation of 
5 Signal energy in a physical system, the attenuated energy is 
usually absorbed by devices such as dampers in laechanical 
systems and resistors in electronic circuits and trans fonr.ed 
to other energy forms such as thermal energy. This may lead 
to undesirable effects and measures, such as using radiating 
10 devices, somstimas have to taken to compensate for these 
effects. When a nonlinear system is employed, instead of 
attenuating che signal energy directly as in the linear 
case, the signal energy at a frequency of interest can be 
spread over a wider frequency band and attenuated by means 
15 of the counteractions between different terms which compose 
the output spectruEi. This means that^ to a certain extent, a 
nonlinear design for signal energy attenuation can reduce 
the requirements for using energy absorption devices and is 
of great benefit in practical applications. 



20 



25 



30 



.n 



Another important application would be, for example, in 
the design of the foundations or the modification of the 
characteristics o£ buildings and structures which are 
earthquake zones. The objective in such an application would 
be to design materials or actuators, tuned as required for 
each structure, which transfer the damaging input energy 
from an earthquake to another more acceptable frequency band 
or spread the energy to be within an acceptable bound over a 
desired frequency range. Spreading the energy using the 
present design should produce significant reductions xn 
earthquake damage. 

(1) Design a nonlinear system to attenuate the energy of the 
signal 
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10 



IS 



20 



25 



U{t) = 



2M., 



;r(a-b) f 



2cos 



It 



t -cosa t -cosJD t 



where d=^3*3, b=^l, M^ = 1.6. 
(2) The spectruja of the signal is evaluated. This is 
implemented by sampling the signal with sampling 
interval T = 0.01 sec to produce 



30 



u(k) = 



2M„ 



k = -199,9r"/0/-v 2000 
and perform a Fast Fourier Transform (FFT) ' on this 
series. The result of the FFT is shown in figure 25, 

(3) Evaluate the range [a.b] of frequencies in the signal 
based on the compured spectrum. The evaluation gives 

a^l,0996 b==3.1416 

Jhis is because the computed spectrum indicates 

lu( jw)|sjO. OS for w «[!. 0996, 3,14X6] 

and 

ja(jw)j<O.OS for other w 

(4) Assume Tihat, as part of the design, output frequency 

range is [c,d] = [o, 7.3] and the required bound over this 

frequency band is Y^'(w) = l. 

(5) Following the same steps for the design of the sysrem 
structure and parameters as in example 1 above, 
Nq, N^and C^, n-N(j,"',N were determined as follows. 

(i)N, 

Nq is determined to be = 1 . This is because 
t.^ =[a,b] = [l.0996, 3.14l6]e[0, 7.3]-[c, d] 
where f^^ denotes the output frequency band contributed 
by the system linear part^ and part of the selected 
output frequency range falls into the linear output 
frequency range t^^ . 
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(iii N 

In this case, it can be obtained using a=l,0990 and 
b=3.1416 that 

f^^ ==[2,1992, 6,2832][J[0; 2.0420] 
f,^==[0, 9.4248] 

Thus, if the maximum noaiinsar order is taken to be 2/ 
then the output frequency range of the system is 

If the maxiinum order is taken to be 3, then the outpui: 
frequency range 

f,U=f,^U^^=[0' 9.4248] 
. Clearly / - f^\^„^ contains the whole selected output 

frequency range [c,d} = [o, 7 . 3] , n is therefore determined 
to be M=3. 

(lii) C^, n = l,2,3. 

C^, Ii = lf2f3. are determined to minimise the 
following expression 



under the constraints at C^>Q, 1 = 1,2/3. 

Notice that for this specific example, 
M = 4000 



i^ = round 



cMT 



27r 



^ round 



0x4000x0*01^ 



fdMT' 



7.6>< 4000x0.01 



2fr 



= 48 



■ round 

Xhe solution to this minimisation problem is 
^C^^O, C2=:2.1932, C3=6.0550 
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(6) Select the NARX model parameters 

based on the results obtained m (5) as 

c:^(l,l,l)=3.Q550, q,3(2,2, 2)=-3. arjd other are zero 

and consequently construct a NARX model as 



20 



25 



- it c,,(i,,gnu(k-i.> ill c„,(i,i,,gflu(k-i,) 

:=c^(l,l)u'(k-l)+q,(2,2)unk-2)+Q,,(lAA)u^k-l)+c,3(2,2,2)u^k-2^ 
=ri,i932u^k-l)-u'{ic-2) + 3.0550u^{k"l)-3u^k-2) 

Notice that the selected NARX model parameters 
15 satisfy the relationship 

i|c«.(i.|=o=c, 

1,-1 

and the different signs selected for 
Gg^d,!) and 00^(2,2} 

and for 

c^3{l,l,l) and c„(2,2,2) 
are such as to gxve effect to the intra-kemel and 
inter-kernel interferences, which is to attenuate the 
energy of the input signal. 
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The frequency domain response of the constructed model 
to the sampled series of the input signal is shown in figure 
27. It can be seen T:hat excellent energy attenuation has 
vjeej^ achieved by the designed system. It will be 
appreciated -hat the xnput energy in figure 2 6 has been 
spread over the designed frequency b5.nd by the nonlinear 
filter. 



10 



13 



Example 3 

This example shows another application of the present: 
irsvention to: the attenuation of signal energy over unwanted 
frequency bands using designed nonlinear ef f ecrs . The 
example also illustrates the effect of the same design on a 
different input signal to dercionstrare the effectiveness m 
energy attsnuarion of the deaigned syateia in dlfferen^c 
Gxrcumstances , 

(1) Design a nonlinear system to attenuate the energy of the 
signal 

1 sin a t —sin fi t 
20 u(t):.2M.-- ^ 

with a^3,3, >ff = l. and M,=1.S. 

(2) The spectrum of the signal is evaluated by sampling T:he 

signal with sampling interval T = 0.01 sec to produce 

1 sin dcT-sin /? 3cT 1 son 3.3xQ.Qlk-s in Q.Ql k 
^^^^^--^ S =2.1.6.-. 

25 ^ ]c = -'1999r"r 0,---/2000 

and then performing a Fast Fourier Transform (FFT) on 
the obtained time series. TUe result of tne FFT is T^ue 
same as that shown in figure 16. 
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(3) Evaluating tue range [a,b] of frequencies in the signal 
based on the computed spectrum gives 

a:=0.6283 b^3,7699 

a.s the coniputed spectrum indicates 
5 lU{jw)1^0.05 for w-e[0,62S3, 3.7699] 

and 

ju(jwi|<0,05 for other w. 

(4) Assise that, as jiiart' of the design, it is desired that 
10 the ou:;put frequency range is [c^d]^[0, 10 .3] and the 

required bound over this frequency range is y^*(w)=l, 

(5) roilowing rhe same steps for the design of the system 
structure and parameters as xn example l above, the 
parameters n^, N, and C^^, n^ljQ/'-vt^ are determined as 

15 follows . 

(i) 

Ng is deteCTined to be Nq=1 due to 

f,^ = [a, b] = [0.6283, 3,7699]e[0, lQ,3] = [c, d] , 
20 which indicates that part of the selected output 

frequency range falls into the linear output 

frequency range f^^ . 

(ii) N 

In this exaniple, N is obtained using a=0.6283 and 
25 b«3.7699. therefore, 

t,^=[0, 7.5398] 

f,^^[0, 11,3097]. 

Hence, if the maximiair. nonlinear order is taken to be 
30 2, the output frequency range of the system is 

^.in=.-f.,U4-{0. 7,5398], 

If the maximum order is taken to be 3^ th© output 
frequency range . 

35 ^|n-J=^,U^.=[0' 11.3097]. 
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As fjjj^,3 contains the whole selected output frequency 

range [cyd] = [o, 7 , 6] , . for the design, is determined 
to be N=3. 

(ill) C,, n = l,2. 3 

C^t n = l,2,3 are determined in such a manner as 
to minimise 



under the constraints of C^^O, i = l,2,3. In this case 
aiao M-4000 but 



ig = round 
i . =a round 



gMT 



= round! 
round 



L 2r I 
10.3x4000x0.01 



= 66 



. 2n i 
" dMT l 
2n- J 

The solution to this minimisation probiam Is 
C,^0, C.=0»2929, 0^=0, 9763 
(6) The NAKX model parameters 

Kf^o^i^u'-.l,) r I,-lr-,K,,-,l,-lr-,K,, for n = N,,^..,N, 
are selected based on the results obtained in (5) as 
= K2 = K3 = 2 

c^(l,l>=0.1923, q,j(l,2)=-0.1, c^l2,l)=0, c^{2,2)=Q, 
<^3(l,Xrl)=0-67S3, qjj(l,2,2)=-0.3, and other q,j(-)'s ars zero, 
and consequently a HJiRX model is constructed 

VPiiXQ li=L L,=l L»l 

• Co2a,l)uNk-l)+CQ^(l,2)u(k-l)u(k-^2)+c^3(l,l.l)u'(l<:-l) 
-hc^j(l,2,2)u(k-l)u^k-2) 

-0.1928u'(k^l)-0.1u(k^l)u(k-2)-hO,CT63u'(k-i)-0.3a(k-l)u'(}c-2) 
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The selected blARX modal paramerers satisfy the 
reiarionships 

iK,(iO|=o=c, 

lIiq«Mhk(-'l|+k<l'2| =!0.192^+|-0.^ =0.2928=q 

and the different signs selected for 

Cog{lyl) and Co2(i,2) 
and for ' - 

c^3(l,l;l) and c„(l,2,2) 
10 are in order to give effect to the intra-kernel and 

inter-rkernel interferences, which is to attenuate rhe 

energy of the input signal. 

The frequency response or the constructed model to the 
15 input signal specified above is shown in figure 28 which 
indicates that the required energy attenuation has been 
realised. 

Thft frequency response of the design above to the input 
20 signal in example 2 gives tha result shown in figure 29, 
The nonlinear system design ^bove clearly also works for the 
signal in example 2 in energy attenuation although the model 
was not specially designed for this signal. This is 
reasonable since rhe magnitude of the spectrum of the signal 
23 in example 2 is less than the magnitude of the spectrum of 
the signal in this example over almost ail of the input 
frequency ^band and over the other f requancy bands the 
magnitudes of the spectra of the two signals are ail zero. 
This illustrates that the above design is effective not only 
30 for the input based on which the design is implemented but 
also for other inputs with magnitude freq^aency 
characreristics less than the magnitude of the spectrum of 
the considered input. 
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XXI. FLOWCHARTS 

Referring to figure 30^ there is shown a flow chart 
3000 depicting the steps of an embodiment of the present 
5 invention which rhe output spectrum is specified over a 
given range of output frequencies , 

The signal to be processed and the desired frequency 
response of the non-linear system to designed are input 
10 via steps 3002 to 300^* ' 

At- step' 3002/ a digital input signal {u(k:) } and its 
sampling' interval T are to be provided. The range of output 
frequencies [c^d] over which the energy of the input signal 
15 is to be transformed is given in step 3004. The output 
frequency range is specified using beginning and end 
frequencies c and d respectively. The distribution of the 
«nergy over the output frequency range [cyd] is re<2Uested or 
specified in step 3006- 

20 

The frequency characteristics of the input signal are 
determined at steps 3008 and 3010. More particularly, the 
frequency components of the digitised input signal/ {u (k) ( / 
are calculated using a Fast Fourier Transform at step 3008. 
25 The range of frequency components contained within th© input 
signal* is determined from the FFT at step 3010, 

Referring to steps 3012 to 3024, the orders of the 

nonlineariries required to realise a desired nonlinear 

30 system and, hence, the energy transf orznation, are 
calculated. 
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Variables n and are set to one and sero respectively 
at 5tep 3012. Thft output frequency components resulting from 
the nth-ordar of system nonlinearities are determined at 
steps 3014 and 3016. 

s 

A determination is raade ar srep 301B as to whether or 
not Nq=0, that is to say^ whether or not the smallest order 

of systesi nonlinearities which makes a contribution to the 
desired energy transformation has been determined. If la 

10 equal to zero, a determination is made at step 3020 as to 
whether or not part or the specified output firoquency range 
falls within fyn* determination is negative, 

processing continues ac step 302S a.t which the value of n is 
increased by one. However, if the determination is 

15 positive, the value of Nq is set to equal n at step 3022 and 
processing continues at step 3026. 

If is not equal to ^ero, then a determination is 
made at step 3024 as to whether or not the specitied output 
20 frequency range lies completely within fv = fvji ^ fYn-l- ■^^ 
the determination is positive, processing continues with 
step 3027. However, if the determination is negative, the 
value of n is incremented by one at step 3026 and processing 
continuas at step 3014. 

25 

Steps 3028 to 3034 detennine the values of the lags of 
the nonlinear model for all values of n = ^q, isfQ+i/^-./N ar.d 

determine the parameters of the nonlinear system to be 
designed, 

30 Step 3036 represents a step for fine tuning of the 

output or frequency response of the designed non-linear 
system. 
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Steps 3038 and 304C deteraine whether or not the 
resulting frequency response is sufficient to meet the 
design requirements. If the frequency response of the 
designed filter is acceptable, the discrete version of the 
5 filter is output at step 3042, Siowever, if the assigned 
non-linear system does not have an acceptable frequency 
response/- then st«p 3044 increments the value of by one 
for all n= Nq, ^Tq+i, ...,N and steps 3030 to 3040 are 
itfera-cively repeated. 

10 ' - ' 

The determination in step 3040 as to whether or not the 
de"signed filter is acceptable is made, for example, by 
calculating the difference between the required output 
spec::rura and the real output spectrum of the designed systera 

15 over the frequency band [c,d] . If the modulus or tne 
difference is below a predeteircainabie threshold over all the 
output frequency band [c,d], then the designed nonlinear 
system may be deemed to be acceptable. However, if the 
modulus of the difference is greater than the threshold at 

20 any frequency over [c,d], then the design is refined. 

Referring to figure 31/ there is shown a flow chart 
3100 for implementing computer code according to a second 
embodiiaent q£ thft presttnt: invention* 

25 

The signal to be processed and the specification for 
the desired bound t:o be imposed upon the output: signal 
frequency charactaristics is input at steps 3102 ro 3106. 

30 At step 3102, a digitised input signal and its sampling 

interval T are to be provided. The rang© of outpur 
frequencies [c,d] over which the energy of the input signal 
is to be transformed is specified in step 3104. A bound 
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Y^* (w) for th« distribution of the energy over the output 
frequency band [c,d] is input at step 3106. 

The frequency characteristics of the input signal are 
5 cLetemitiftd at step 3103 and 3110. More particularly; the 
frequency components of the digitised Input signal^ {u (k) } ^ 
are calculated using a Fast Fourier Transtorm at step 3108. 
The range of frequency components contained within the input 
signal i$ determined from the FFT at step 3110. 

10 

Referring to steps 3112 to 3124, the orders of the 
nonlinearities required to realise a desired nonlinear 
system and/ hence; th6 energy transformation, are 
calculated, 

15 

* Variables n and Nq are set to one and zero respectively 
at step 3112. The frequency components resulting from the 
nth-order of system nonlinearities are determined at steps 
3114 and 3116. 

20 

A dfttemiination is made at Step 3118 as to whether or 
not Up = 0, that is to say, whether or not the smallest 
order of system nonlinear it ies which contributes to the 
desired energy transformation has been determined. If is 

25 equal TO zero, a determination is made at step 3120 as to 
whether or not part of the specified output frequency range 
falls within fyn* ^If the determination is negative, 
processing continues at step 3125 in which the value of n is 
increased by one. However, if the determination is 

30 positive, the val^je of iJ^ is set to equal n at step 3122 and 

processing continues at step 3126. 

60 



SUBSTITUTE SHEET (RULE 26) 



30-08-00 12:05 From: HARRISON GOODARD FOOTE 



^41 142730312 



T-663 P 68/99 Job-687 



WO 99/45644 PCT/GB99/00550 

If, at step 3118, it is detecmined uhac Nq is not equal 
to zero, then a determination is made ar step 3124 as to 
whether or not the specified output frequency range lies 

corttpletely within fy = fyn ^ ^Yn-l* detemination is 

5 positive, processing continues with step 3127. However, if 
the determination is negative, the value of n is incremented 
by one at step 3126 and processing continues at step 3114, 

Steps 3123 to 3134 determine the values Cj^ for all 
10 values of n = Nq^ Nq^^^/.^wN/ which represents- the summation 
of, the rTLOdulus of the values of "he parameters of the 
nonlinear ' system. 

At step 3136 Che paramerers of a NARX model are 
selected given the const:ra.ints imposed upon the non-linear 

IS system obtained in steps 3128 to 3134. A determination is 
made via steps 3138 and 3140 as to whether or not the 
frequency response beyond the frequency band [c,d] of the 
designed filter is acceptable. If the frequency response is 
acceptable, the filter design parameters are ourpui: ai: step 

20 3142. However, if the filter charactsristics are not 
acceptable at frequencies outside the frequency range [c,d]^ 
then a conventional filter, H(q"^), is designed, at step 
3144, in order to reduce or obviate the frequency response 
of the designed filter outside the range of frequencies 

25 [c,d].^ Finally, the design is completed at step 3146 by 
combining the designed non-linear and linear filters, if 
any* 

The detenoination m step 3140 as to whether or not the 
30 frequency response of the designed filter is acceptable 
outside the range of frequencies [c,d] is made, for example, 
by comparing the modulus of the frequency response beyond 
the frequency range [c,d] with a p rede terminable threshold. 
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If the modulus is below the threshold over all the frequency 
range outside [c,d], then the designed nonlinear system may 
be deemed to be acceptable. However^ if the modulus is 
greater than the threshold at any frequency beyond [c,d], 
5 the design is refined, ■ 
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XV, THME MORE DESIGNS EXftMPIiBg 

IV g DESIQ? OF NONLIKSi^ FILTERS WITH SBECIFICATIONS FOR 
BOTH THE MRGtrrTtrDH AlID PHASE OF OtTTgUT FREQUENCY RESPONSES 

5 

It will be appreciated that the basic principles o£ the 
present invention can be applied to the design of nonlinear 
filters based on specifications for both nhe magnitude and 
phase of output frequency responses . The ability to 
10 modulate phase as well as or instead of magnitude is of 
particular importance in teiecoinmunicacion applications. 

Consider, a filtering problem such that given an input 
spectrum U{jw) over a frequency band [a,b] and a desired 
output specrrimi Y\jy) over anocher frequency band [C/d], a 
15 filter is required zc be designed so that the output 
frequency response y{j^) can match the desired spectrum 

^'(j^) closely as possible in terms of both magnitude and 
phase characteristics. The basic principles in Section I can 
be directly applied to realise the design of such a filter. 
20 The procedure to be followed is described below. 

First, a nonlinear filter 

is designed using the basic principles descrii^ed in Steps 
(i) — <iv} and (v.l) Part 1 in Section I to produce a 
25 frequency response Y^ijw) such that Y^Uw) can match Y^\j^) 
over the specified output frequency range [c,d] as closely 
as possible in terms of both the magnitude and phase. 

Secondly, if required, design a linear filter with a 
frequency response function H^{jw) such that, ideally. 
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Tbxs can be used to ojnprove Y^(jij) , the result obtained from 
ths nonlinear design, $o that the frequency response of 
the corresponding output 

provides a better match to the desired spectrum V'ij^) * 

Thirdly, if required, design a linear phase FIR (Finite 
Impulse Response) band pass filter H^U^) ^^t:h the ideal 
magnitude frequency characteristic 

[O otherwise 

10 and linear phase over the frequency range. Then construct 
the designed filter using linear filters H,(jw) and H,U^) 
and nonlinear filter N[uit)] as shown in figure 32 zo yield 
the output frequency response 

15 The second and third steps above follow the design 

principle described in Step (v.l) Part 2 in Section I so as 

CO augment the performance of the nonlinear filter 
yt(fc)=== N[u(t)] designed in the first step. 

Ideally, ■ yCj^^)-H,{jp^)r,(jv) implies that 

} otherwise 



20 1^0^)1== 
and 



due to the linear phase characteristic of H.^D^) where ic, is 
a coefficient which is a function of the order of Hj(jw) . 
25 This indicates that the output frequency response of the 
designed filter is ideally the same as the desired response 
over the output frequency range [c,d] except for a linear 
phase difference between the real and the desire4 phase 
response. This is in fact an unavoidable phenomenon if band 
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pass filtering such as the effect of H^ijw) is applied for 
the design. However, this is still an ideal characteristic 
since a linear phase implies that the filter group delay is 
constant, the filtered signal is simply delayed for some 
5 time determined by k^, and the wave shape of the processed 
signal is preserved, that is, there is no phase distortion 
in the filter output. 

The results of a specific nonlinear filter design 
obtained using the above procedure and based on 
10 specif icaticns for both the magnitude and phas^ are shown m 
figures 33-35. 

The' given input for this design was produced using a 
white noise sequence uniformly distributed in [0,4] and band 
limited within the frequency range [a,b]=[l,S] under a 
15 sampling interval 0,02 s. The magnitude of the spectrum of 

the given input is shown in figure 33, 

•p^e desired output spectrum was chosen to be 
|exp(-500w)^j(600^--) n^.^^i^-j 

y = \ 1000 L i L -> 

[O othexwise 
Figure 33 also shows the comparison between the 
20 magnitude of the output spectrum Yijv/) of the designed 
nonlinear filter and the magnitude of the desired spectrum 

Figure 34 shows the comparison between the phase angle 
of r3(jv),"1:he output spectrum before linear phase filtering, 
25 and the phase angle of the desired spectrum Y\jw) . 

Figure 35 shows the phase angle of the applied linear 
phase filter . 

It can be obseinred £roiii the above that a nonlinear 
filter designed using the basic principles of the present 
30 invention can produce an output frequency response which 
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Figure 35 shows the phase angle of the applied 
linear phase filter . 

It can be observed from the above that a nonlinear 
filter designed using zhe basic principles of the present 
5 invention can produce an output frequency response vhich 
satisfies a design specification in terms of both the 
magnitude and phase. 

The discrete time model description of the designed 
filter is given below where the nonlinear part of the 
10 model is the discrete time model description of the 
nonlinear filter 

■y,(t)-NWt)] 

and the linear parr of the model is the discrere time 
model description of the linear filcer, tha fre<:jaency 
15 response function of which is 

Nonlinear part of the raodel: 

20 y(k)= +(3.123e+06) u(k-l) +(-2,2446+07) u(k-2) 4* ( 6 . €02e+07) 
u(k-3> 

+(-l.027e+08} u(k-4) +{8.961e+a7) u(k-5) +(-4.l8e+07) 

u(J<-6} 

-t'(8.l76e-*-06) u(k-7) +{4,235e+Q9) uCk-l}u(k-l) 
25 +{-l,772e-J-10) u(k-l)u{k-2) +(5.a97e+05) u(k-l}u(k-3) 

*-H{7.213e+09) u (k-l) u +(-3.66e+09) u(k-l)u(k-5) 
+ L?1.869e-f08) u(k-l)u(k-6) + (6,736e+-07) u (k-1) u (k-7). 
+ (--6.777e+08) u (k-2) u ( k-'2) +(8,003enO) u(k-2)u(k-3) 
+ (-7.57e+10) u(k-2)u(k-*4) +(6.023e+09) u(k-2)u(k-5) 
30 +(a.648e+09) u(k-2)u(k-^) + {-4 . 365e+03 ) u(k-2)u(k"7) 

+ (-8.365e+10) u(k-3)u(k-31 +(2.268e+10) u ( k-3) u ( 3^4 ) 
+ {l,OOSe+ll) u(k-3)u(k-5) + (-3-795e+10) u(k-3)u{k-6) 
+{-2.5166+09} u(k"3)u(k-7) +{1.054e+ll) u(k-4)u(k-4) 
+ {-l,97e+ll) u{k-4>u(k-5) +(l,109e+10) u ( k-4) u (k*-6) 
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r (1.907e4-i0) u(k-4)u(k:-7) T{2.352e+I0) u(k-5)u{k-5) 

+ (8.i73e+lQ) u(k-5)u(k-6) + (-3 , 318e+10) u ( Jc-5) u (k-7) 

+c~4.218e+10) u(k-6)u(k-6) +(2.042e+10) u(k-6)u(k-7) 
+(-1.658e+09) u(k-7)u(k-7) 

5 

Linear parr of the model: 



y(k)^b(l)u(k)+b(2)u(k-l)f, . . , , +b (m) u ( k-m) -a (1) y (k-1) . . 
a (n) y (k-n) 

10 with t].-2 and m=303 where 

[a{l) ,a(n)] - 

-1.8 6215871605398 0,95715116734794 

[b{l}/ b(m)] - 

15 l.Qe-03 * 

Columns 1 through 4 

0. 00000003328202 0.00000021432147 0 . 000QOO80368128 

0,00000227704991 

Columns 5 through 8 

20 0 .0000053787 6 0 . 000011134 9910 S 0.0000207 9938253 

0.00003571014847 

Columns 9 through 12 

0.0000 570 4570668 0,00008 54 3478 679 0.00012079738 459 

0.00016141663710 

25 Columns 13 through 16 

0*00020406372559 0.00024351205626 0.00027258483896 

0,00028246765302 

Columns 17 through 20 

0.0002,6339508245 0.00020572744049 0,00010137844559 
30 -0. 00005450964923 

Columns 21 through 24 

-0. 0002 6181250429 -0,000513 818 00480 -0,00079 60337 3 996 
-0.00108568381434 
Columns 25 through 28 
35 -0.00135214637123 -0.00155851835832 -0.00166438952001 
-0,00162976908378 
Columns 29 through 32 
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-0.00141994747987 -0.00101090816448 -0.00039475153522 
0,00041552374840 

Columns 33 through 3 6 

0. 00138259540242 0.00244344402372 0,00351062196833 

5 0.00447655873267 

Columns 37 through 40 

0,00522101718300 0.00562147075427 0 . 0055 5573253 464 

0.0049S62727fi545 

Columns 41 through 44 

10 0.00377365278712 0.00198963381599 -0.00032378155160 

-0.00303968275897 

Columns 45 through 48 

-0.0059 676904053 9 -0,0088 6230677976 -0,011428433 94 688 
-0.01339358387193 
15 Columns 49 through 52 

-0.01443551408 603 -0 . 01431325012233 -0.01284873916639 

-0.00996641286493 

Columns 53 through 56 

-O.Q0571602453045 -0.00028726506096 0.0059884125 6370 

20 0.01264965179911 

Columns 57 through 60 

0.01913182739200 0 . 0248 0771068003 0.02903991461538 

0.03124135106508 

Columns 61 through 64 

25 0.03093862246516 0.02783233132424 0.0218478 68 63333 

0.01317044094614 

Coluxnns 65 through 68 

0.002258 96524855 -0,01016501602446 -0.023155158 92661 
-0,03560414378179 
30 Columns 6^9 through 72 

-0.0463300789195S -0.05417 343987784 -0.05813822823951 

"0.05744137376473 

Columns 73 through 76 

-0.05166958483783 -0,04082303832663 -0.02536256454818 
35 -0.006Z1129546254 

Columns 77 through 80 

0.01528703636150 0,^03745307958913 0,05839584 942467 

0.07616729599331 
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Coluians 81 through 8 4 

0,oa9 93S37701557 0.09516078483677 0,09376093809767 
0.OB42447S625381 

Columns 85 through 88 

5 0.06630350720636 0.04234527006783 0 . 0124S5€$4280lS 

-0,02064087852648 

Columns 89 through 92 

-0,054346138 62312 -0.08578337461088 -0.11209316760456 
-0.13067820746590 
10 Columns 93 through 96 

-0,13945154456062 -0 . 13705390892471 -0,123018230477 38 

-0.09736316009190 

Columns 97 through 100 

-C. 06310299056015 -0.02116847625143 0.02 475895107438 

IS 0.07098889554710 

Columns 101 thi^Qugh 104 

0.11362133287669 0.148 87 942949652 0.17344622849594 

0*18477606135825 

Coiunins 105 through 108 

20 0.1313498 4027062 0.16284992659307 0,13023456167 4 87 

0.08570350882143 

Coliiitins 109 through 112 

0.Q3255O64587105 -0.02508361142904 -0 . 08254318260252 
-0.13504269472306 
25 Columns X13 through 116 

-0,17807650738416 -0.20781448143133 -0.22144885867665 
-0. 21746018536730 
Coluians 117 through 120 

-0.1957768088465S -Q . 15781155123986 -0.1063 6991902139 
30 -0,045435 6718 5132 

Columns 121 through 124 

0.02014930576277 0.08508414798056 0.144053048 99329 

0.19217566873157 

Columns 125 through 128 

35 0,22542441249903 0.24097171607669 0.23743645368538 

0.21500740460840 

Columns 129 through 132 - 

0 .1754325 6593 954 0 . 1218 7 4992 477 61 0 . 058 647 63 64 003 5 
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-0 .00914923229050 
Columns 133 through 136 

-0,076057259374 66 -0.136716636833 33 -0.19 631209098 366 
-0.220968693S50S9 
5 Colurons 137 through 140 

-0.23S064158^5179 -0.23 6431475X1610 -0-21643603250516 

-0.17992231311682 

CoiuKins 141 through 144 

-0. 13003701658436 -0,07094606546323 -0,00747191949549 
10 0.05531598034943 

Columns 145 through 148 

0.1125214 9812387 0 ■ 15981666793597 0.19377819943088 

0\2121387423''9435 

Columns. 149 through 152 

15 0.2139342 9471004 0,19953378 954962 0.17058700949905 

0.12979660938955 

Columns 153 through 156 

0.0807083258 9097 0.02736960335274 -0.02 600882989955 
-0.07539394440027 
20 Columns 157 through i60 

-0 » 1172 4050322409 -0 . 14 87 4 67 98 8 0297 -0,16803574024 719 

-0.17425382855726 

Columns 161 through 164 

-0.16758132308460 -0,14916007746448 -0.12094944070741 
25 -0.08552673296083 

columns 165 through 16B 

-0.04595223365109 -0 . 00SQ201703538 4 0.03398322634818 
0.06849574758055 

Coiurans ^69 through 172 

30 0,09633411698109 0.11601725927068 0.12563256434628 

0.12804494388799 

Columns 173 through 176 

0.12075316710299 0.10582180562013 0,08475985344594 

0.05937858562038 

35 Columns 177 through 180 

0.03164147930605 0,003518162377 07 -0 . 023147409599€G 

-0.04674674074459 

Columns 181 through 184 



70 

SUBSTITUTE SHEET (RULE 26) 



•OC 12:07 Froni:HARRISON GOODAfiD FOOTE 



H4i 142730312 



T-663 P r8/9S M-BBJ 



WO 99/45644 PCT/GB99/00550 

-0,06598945884170 -0 . 07995651S94^41 -Q . 0S8 12411270222 
-0. 09036059631364 
Col^jinns 185 through 188 

-0.08690017680633 -0.07829917708583 -0.06537295841595 
5 '0.04913952047510 

Colunnis IS 9 through 192 

-0.03073930999045 -0.01137 008 308799 0.00778109613199 
0.02559935810198 

Coiuinns 193 through 196 

10 0.04109626424759 0,05345544314259 0.06206855688305 

0.06656076607150 ' ^ 

Coiiarais 197 through 200 

O' 06680468888031 0.06292254225138 0,05527568217195 

0.04444169619049 

15 Coluims 201 through 204 

0.031179 46545587 0.01638338990942 0,00102875697948 

•0.01388895099522 

Coiuirms 205 through 208 

-0.(3274172418 9660 -0,03870846951445 -0.04707559699096 
20 -0.05203784736309 

Colmufia 209 i::hxough 212 

-0-05335260858275 -0,0510308 4289941 -0.045334 43158680 

-0.03675527219878 

Columiis 213 through 216 

25 -Q. 02597741626758 -0 ,01382495500742 -0,001199592^3130 

0,01098723154804 

' Columns 217 through 220 

0.023^8822958380 6 0,03075319656401 0.03704059636285 

0.04039599306472 

30 Coiuinns'^221 through 224 

0.04070212260137 0.03807460516092 0.03284505620907 

0.02552746332356 

Coluiuns 225 through 228 

0. 01677103005022 0.00730378144712 -0.00212808430644 
35 -0.01081988478115 

Columns 229 through 232 

-0.01816221252030 -0.02363447649585 -0.02708498340864 
-0.02624561802129 
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Columns 233 through 236 

-0.02723060984541 -0,02427029083261 -Q . 01973202360355 

-0,01408149907660 

Columns 237 through 240 

5 -0,00783828203366 -0.00152 977531123 0 . 00435233046S71 

0.00938952420923 

Columns 241 through 244 

0.01326503243370 0.01578162483958 0,0168 6791315533 

0.01657347478835 

10 Columns 245 through 248 

0.01505411203293 ,0^,0125493 4451451 0 . 0093S476560966 

0.00579214028759 

Colu22ins 249 through 252 

0.002ia0O8037598 -0.00119216975076 -0,004085818 00286 
15 -0.00632711324010 

Coluinns 253 through 256 

-0 .00731483560799 -0 , 00352066157236 -0.00848329376098 
-0.00779740923517 
Coluinns 257 through 260 
20 -0.00 65^902X70090 -0 . 00504908820003 -0.00331691201488 
-0.00156482692685 
Coluinns 261 through 2 64 

0.00006463440339 Q . 00145324730714 0.00253783620427 

0.00326220436512 

25 Colmnns 265 through 268 

0.00362556118310 0.00365304002511 0.00339409451670 

0,00291465798950 

Goiucnns 260 through 272 

0.00228895177653 0.00159174429109 0.00039171267597 

30 0.00024^6506918 

Coluttins 273 through 27 5 

-0 .00030123128510 -0.00072387920912 -0 . 0Q100998 4811S8 
-0.001161S8515783 
Columns 277 through 280 
35 -0,00119315515203 -0.00112532046764 -0,00098439594641 
-0,00079761965498 
Columns 281 through 284 

-0.000S9067279083 -0.0003855777567 6 -0.00019936402666 
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-0,00004349568936 
Coiucnns 285 through 288 

0.0O0O76O24S6265 0 . Q0015B01946385 0.0002051369103 6 

0.00022267635728 

5 Columns 28 9 T:hrough 292 

0.00021739726997 0 . Q0019645030203 0 *00016652267489 

0.00013324716319 

Columns 233 through 296 

0.00010088214188 0.00007223559481 0.00004878213424 

10 0. 00003090906954 

Columns 297 through 300 

0.00001822614592 0,00000988083116 0,000004 83 478353 

0,00000207399807 

Col-ums 301 through 303 

15 0.0000007 4209218 0 . OOQC0Q200733Q4 0 . OOOOQ003164S28 

TABLE 1 
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IV, 2 OSSIGN OF NONLCTEAR rXLTSRS TO FOCUS KtlERGY FROM 
DIFFEKEOTT FREOtJENCY 5AKDS TtrgQ A SIITGLE FHSQUEHCY SAIID 

The above einbodimenrs of i:he present: invention 
5 illuSurate how the principles underlying the mventxon can 
be used to transfer energy from one frequency or frequency 
range co another frequency or frequency range. However, it 
will be appreciared that the present invention can equally 
wsil be utilised to focus energy from a given freouency 
10 band or given frequency bands into a single/ and 
preferably, narrower frequency band. 

The problem to be addressed in this design is that 
given an- inpur spectrum which possesses nonzero magnitude 
characteristics over two different frequency bands [a^,b^] 

15 and [a.,i?.]y where d,>i:>^r a filter is required zo be designed 
to focus energy from the two different inpur frequency 
bands into a single output frequency band [c,d]r where 
c> b^, d<a^^ ; with the spectrum of the filter output over nhe 
output: frequency range [c,d] satiisfyxng certain 

20 specifications. The specxf xcations to be satisfied can be 
specified in terms of magnitude^ phase or both. 

The only difference of this design compared to the 
designs in Section I is that the determination of the 
maximum order M of the filter noniinearities should not 
25 only ensure that the specified output frequency band [c,d] 
is covered by the filter's output frequency range, but also 
should ensure that the output energy over [c,d] Is derived 
from the input frequency range [^i)b,] and the input 
frequency range: [a,,^^] * Therefore, the principles for this 
30 design are eicactly the same as those described in Secrion I 
except for step (iii) which should follow the new principle 
above . 
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Consider an example where a nonlinear filter is to be 
designed to focus the energy from two different frequency 
bands [3,4] and [10,11] of an input signal 



U(t) = IM^ 



27t 



5 where b, = 4, 3, K = 1 1, a, = 10, - 0.03 , into a aingle £::equency 

band around frequency £==7 wirh the single ourput frequency 

band as small as possible and the output magnitude 
characteristic at frequency f=7 being 3M(^=0.09* 

It can be shown that when subject to the excitiation of 
10 an"~ input with the frequency spectrum over nwo different 
frequency bands [a,,bj and [a,,^,]' output frequency 

ranges of a nonlinear system contributed by the system's 
2nd order nonlinearity is 

[2a,^], [23,^], [a, +a„i>, +h.], [O,^, -a,]u[0,i:i, -a.] and [a,-b^,}x -a,] 

15 Su22srirut:lng iij = 4, a, = 3, jb. = a, = 10 into - -bj and b, - a^ 
yields 

[a, -b,,b,-a,] = [6,8] 
Clearly this output frequency range covers frequency f-1 
and the system output energy over this frequency band is 
20 derived from the input frequency range [a^bj = [3,4] and the 
input frequency range [a^^b^] = [lO,ll] . Therefore, the maximum 
order of system nonlinearities can be raken to be N=2 for 
this specific design. 

Moreover, following the same design principles as 
25 described in Section I, a discrete time nonlinear filter 
under the sampling period r^^l/50s can be obtained for this 
specific design* Figure 36 illustrates schematically such 
a filter. 
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Figures 37 and 38 show the input and output frequency 
spectra of this filter respectively and clearly indicate 
that the input energy over the two different input 
frequency bands [3,4] ^nd [10,11] have been, as required, 
5 successfully focused into a very small band around 
frequency f =7 , 

IV. 3 DESIGtl OF SPATIAL DQMajN NQNLXKEAll FILTBRS 

Although the above embodiments have been described 
10 with reference to tne^ design of non-linear time-domain 
systems, it will be appreciated that the present invention 
is equally applicable to non-linear spatial-domam systems. 

In -the one dimensional case, the present invention can 
be directly applied and the only difference of the one- 

15 dimensional spatial domain case over the above time-domain 
embodiments is that the time-domain variables t and k in 
the . continuous and discrete time filter equations become 
the continuous and discrete spatial variables X and in 
the spatial domain filter equations. 

20 Consider a one dimensional spatial domain filtering 

problem that is the same as the second problem addressed in 
Section II. 1 except that the time variables t and k become 
the spatial variables X and in the present design. 

Figure 39 shows the block diagram of the spatial 
25 domaia nonlinear filter designed in this case. It can be 
seen that the spatial -domain filter is exactly the same as 
that shown in figure 9 which is the block diagram of the 
corresponding time domain nonlinear filter designed in 

Section II . 1 except that the time variables t, k, T, ana 

30 s (second) in figure 9 have been changed to the spatial 
variables X, Xj,, AX, and m (meter) in figure 39. 

Figures 40 and 41 show the effects of spatial domain 
nonlinear filter. Again, except for the time variables 
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being changed zo spatial variables, rhe figures are the 
same as figures 11 and 13 which are the results of the 
corresponding time doiriam nonlinear filter. 

The filtering effect m the time domain af the second 
5 nonlinear filter designed in Section II. 1 indicates an 
energy rransfer £rom the input frequencies t^^; = 1 and 5^^2=2 
to the output; frequency w-_jo=4. Figure 4 0 indicates the same 
effect but in the spatial domain. In this case, the spatial 
domain nonlinear filrer transfers energy from, the input 
iO spatial frequencies and w^. = 2 to the ourpur spatial 

frequency 

^When this one dimensional spatial domain nonlinear 
filter is applied for one dimensional image processing, 
where figure 41 represents the intensities over the spatial 

15 variable X of the images before and after the processing, 
the energy transformation effect of the spatial domain 
nonlinear filter can be furrher illustrated by figures 42 
and 43 ♦ Figux^ 42 shows the one dimensional input xmage 
with energy located in spatial frequencies f^^i = 1 and w^2 ~ 2 

20 and figure 4 3 shows the one dimensional output i:mage with 
energy located in another spatial frequency due to 

the effect of nonlinear filtering. 

Although the present invention has been described 
25 above with reference to a one dimensional case in the 
spatial do;main/ it will be appreciated that the present 
invention is not limited thereto. The present invention can 
equ^aily w«ll be applied to the design and implementation of 
non-linear systems for realising transfer of energy from 
30 first m-dimensional spatial domain frequencies to specified 
second n~dimensianal spatial domain frequencies, where m 
and n are greater than one, 
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It will be appreciated that an application cf the 
transfer of eaergy from first m-dimensional spatial domain 
frequencies to specified second n-dimensional spatial 
domain frequencies would be digital image processing or 

5 filtering. In such a case m=n=2 . The present invention 
can be designed to produce filters that ap^rate upon, 
digxt^l iiuages. The filters can be dasigned to perform 
nuznerous different functions such aSx for example, image 
compression or filtering to remove ncise or vary the colour 

10 space of an image. 

The phrase "transfer of energy" includes , without 
limitation/ the processing of a first signal^ specified in 
the time domain or spatial domain> to produce second signal 
having predetermmabie characteristics including a 
15 specified energy distribution. 

Furthermore^ it will be appreciated that zhe terms 
"frequency range" and "range cf frequencies" includes a 
group of frequencies or frequency group. A group of 
frequencies comprises a plurality of frequencies spatially 
20 distributed within an n-dimensional space or over a 
subspace . 

It will be appreciated- by those skilled in the art 
that the phrase '^specified spectrum'" relates to the 
specification of at least one of either th^ magnitude and 
25 phase of a signal or frequency components of a signal and 
IS ecfually applicable to situations in which only the 
rnagnitudS' is specified and to situations In which both the 
magnitude and phase are specified* 

It will be appreciated that the input and output 
30 signals of any or all of the various embodiments of the 
present invention may be time or spatial domain continuous 
or discrete input or output signals. 
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Referring to figure 44 there is shown a data processing 
system 4400 upon which einbodiments of the present invention, 
that is to say/ rhe non^linear systeras and the methods of 
design thereof can be implemented or realised. It will also 

5 be appreciated that zh^ non-linear systems so designed may be 
implemented in digital form using the data processing system 
or other suitable hardware and software. The data processing 
system 4400 comprises a central processing unit 4402 for 
processing computer instructions for implementing the design 

10 of nonlinear systems , given an input signal and particular 
outpur signal requirements. The data process'ing system also 
comprises a msmciry 4404 for storing data to be processed or 
the results 'of processing as well as computer program 
instructions for processing such data, a system bus 4406, an 

15 input device 4408, an output device 4410 and a mass storage 
device, for example, a hard disc drive 4412. 
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1, A inethod for designing a non-linear system for 
transferring energy from a time or spatial doiaam 

5 input signal having a first spectrum at a first pre- 

determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-^de terminable frequency or 
range of frequencies, said merhod comprising the 

10 steps of 

identifying or specifying the first spectruea o£ the 
time or spatial domain input signal from which energy 
is to be transferred/ 

specifying the second spectriim of the time or spatial 
15 domain output signal to which said energy is to be 

transferred, and 

calculating/ usiag a frequency domain descriptxon of 
said output signal, for example, the output spectrum, 
expressed m terms of a frequency domain description 

20 of sAid input signal and coefficient: s qz s ^ime or 

spatial domain description of a generalised non- 
linear system, the coefficxents of the time or 
spatial domain description of the generalised non- 
linear system in order to give effect to the energy 

25 transfer, 

2. A method as claimed in claim 1, further comprising the 
step of 

selecting a time or spatial domain description of the 
30 generalised non-linear system; 

determining or defining a frequency domain description 
of the time or spatial domain input for the generalised 
non- linear system; ^nd 
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determining or defining zhe frequency domain 
description of the output signal, for example, the 
output spectrum, of the generalised non-linear system 
expressed m terms of the frequency domain description 
of said input signal and the coefficients of the time 
or spatial doraaxn description of a generalised non- 
linear sysT^ein* 



3. A method as claimed in either of claims 1 or 2, wherein 
10 the frequency domain description of t:he input signal is 

U(jw)i the time or spatial domain description of said 
generalised non-linear system is given by zhe 
generalised NARX model 



15 where 

yjk)=i: t c„(i,,...,i..,)ny(ic-i,)ffu(k-i,) 

with 

5f„ K„ 

p + q = n, = ..;K^, i = l,...,p + q, and £ - S"' S 



{C2) 



20 



the frequency domain description of the output of the 
generalised non-linear system is given by 

where 



25 N is the maxiraum order of dominant system nonlinear it ies, 
{(Odcr^ 



81 



SUBSTITUTE SHEET (RULE 26) 



30-08-00 12:10 FromiHARRlSON GODDARO FOOTE Hi]UZr3Qm T-6S3 P BBm Job-687 



WO 99/45644 PCT/GB99/00550 

deuotes an integration over the nth^dimensionai hyper- 
piane w,t,'",+w^ = w , and jw,,"vjw^) , n = I,...,I are 
generalised frequency response functions of the non- 
5 linear system. 

4. A method as claimed in any preceding claim, further 
comprising the step of 

determining a mapping between the time or spatial 
10 domain description of the generalised nonlinear system 

and the- frequency domain description of the generalised 
nonlinear system, 

5, K method as claimed in claim 4, wherein the mapping from 
15 the time or spatial domain description of i:he 

' generalised non-linear system to the frequency domain 
description of the system is given as 

|l-£c,,(li)exp[-j {w,+,-,+w^) lJ|H^(jw,,"vjwJ 



ZO (C5) 
where 



(C6) 
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A method as claimed ir. any preceUing claim, further 
comprising the steps of 

defining or determining a general relationship between 
the input and ourput frequency or frequency ranges of 
the generalised non-linear system. 



10 



15 



20 



25 



7 » A method as claimed m any preceding claim, wherein the 
relationship between the input and output frequencies 
or frequency ranges is given by the following 

"%=£f.U£.,., (CI) 
where denotes the range of frequencies of the outpun/ and 
f.^ and f^^^^ denote the ranges of frequencies produced by the 
Nth-order and ( N --1) th-order nonlinearities, and 

nb 



>7hen 
when 



(a + b) 
nb 

(a + b)" 



na 
(a + b) 
na 



<1 



>1 



ICS) 



(a + b) 
n=N and N-1 
where I,] relates to or means take the integer part, 
na 

1 = +1 

L U-hb) J 

Ij, =[na-k{a-h-b) ,nb-k(a+b)] fox k = 0,...,i' - X, 
l^. =[0,ub-i'{a+b) 



and the frequencies of the signal to be processed are 
in the range [a,b] and given [a,b] and the required 
output frequency range fy, the method further comprises 

the step of determining the smallest N from the 
relationship above for the generalised non-linear 



83 

SUBSTITUTE SHEET (RULE 26) 



30-08-00 12:10 FrorHARRfSO^f G0DOARD FOOTE 



+441142730312 



T-563 P Sl/99 Jab-687 



WO 99/45644 PCr/GB99^00S50 

systein vhich. can bring about the specified frequency 
domain energy transformation. 



8. A method as claimed in claim 7, wnerein, with M having 
5 been detenainsd and K^,n= 1,,..,N , bea-ng given a priori, 

the raethcd further comprises the steps of; 

taking N as N and determining the coefficients of the 
time or spatial domain model of the generalised non- 
linear system from th© description for the system 
10 output spectrum given in terms of the spectrum of the 

input signal and the coefficients of the time or 
spatial domain model of the said generalised non-linear 
system* 



15 9, A method as claimed in claim 8, further comprising the 
steps of 

substituting j^i/'-'/Dw^) given in (C5) into (C4) , and 
substituting the resultant expression for YJjw) into 
(C3) to obtain the description for the system output 
20 spectrum in terms of a function of the spectrum of the 

input signal and the coefficients of the time or 
spatial domain model of the said generalised nonlinear 
system. 

A method for realising or manufacturing a. non^-linear 
ayarerci for transferring energy from a time or spatial 
domain input signal having a first spectrum at a first 
pre-determinable frequency or range of frequencies to a 
time or spatial domain output signal having a second 
spectrum at a second pre-determinable frequency or 
range of frequencies/ the method comprising the steps 
of 



25 10. 



30 
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(a) designing the ncn-iinear system using the method as 
claimeci in any of claims 1 to 9; and 

(b) materially producing the non-linear sys-cem so 
designed or using the non-linear sys-cem so designed to 

5 modify materially the transfer function of an existing 

linear or non-linear system. 



11. A data processing system for designing a non-linear 
system for transferring energy from a rime or spatial 
10 domain input signal having a first spectrum at a first 

pre-determinable frequency or range of frequenciea to « 
time or spatial domain output signal having a second 
spectrum ac a second pre-determinable frequency or 
range of frequencies, said sysx:em comprising 

15 means for idenrifying or specifying the first spectrum 

, of the time or spatial domain inpur signal from which 
energy is to be transferred, 

means for specifying the second spectrum of the time or 
spatial domain ourput signal to vhich said energy is to 
20 be transferred, and 

means- for caicularing/ using a frequency domain 
description of said output signal, for example, the 
output spectrum, expressed in terms of a frequency 
domain description of said input and coefficients of a 
25 time or spatial domain description of a generalised 

non-linear system, the coefficients of the time or 
spatial domain description of said generalised non- 
linear system in order to give effect to the energy 
transfer. 



30 



12. A system as claimed in claim 11, further comprising 
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means for selecting a time or spatial domain 
description of the geaeraiised non-linear system; 

means for determining or defining a frequency doiaain 
description of the time or spatial domain input for the 
5 generalised non-linear system; and 

means for determining or defining the frequency domain 
description of the output of the generalised non-linear 
system expressed in terms of zhe frequency domain 
descriprion of said input: signal and the coefficients 
10 of the rime or spatial domain descripri-on of a 

generalised r.on- linear system. 



13, A system as claimed in either of claims ii or IZ, 
wherein the frequency domain description of t:h« input 
signal is V(jw}t the time or spatial domain description 
of rhe generalised non-linear system is given by the 
generalised NARX model 



15 



where 



20 y„(k)=z i c,,(i,,...,i,,,)ny(k-ijnu(ic-i,) (CIO) 

with 

p + q = n, 1, K^, i = l,..,,p4-q, and £ 



the frequency domain description of the outpur of rhe 
25 generalised non-linear syatem is given by 

Y(jw) = SY,(jw) (Cll) 
where 
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^^^^ «1+^.,>«^:*W 4=1 

N is the ma^imuxa order of dominant system nonlinearities; 



5 denotes an integration over the nth-dimensional hyper-piane 
w^-f.-'.+Wjj - w , and j w^^, jw^) , n-l,.*.,^J are generalised 
frequency response functions of the said non- linear systeni. 



14. A sysrem as claiiaed in any of claims 11 to 13, further 
10 comprising 

means for determining a mapping between the tim^ or 
sparial domain description of the generalised nonlinear 
^system and the frequency domain description of the 
generalised nonlinear system. 



13 



20 



15, A system as claimed in ciairci 14, wherein the mapping 
from the time or spatial domain description of the 
qeneraiised non-linear system to the frequency domain 
description of rhe system is given as 



= 2l<^^.a^-Aj exp[- j(wil,+,-,+w„lJ] 



, (C13} 

where 
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(C14) 



IS. A system as claimed in any of claims 11 to 15, further 
3 comprising 

means for defining or detenraning a general 
relationship between the input and output frequency or 
frequency ranges ^of the generalised non- linear sysTiein. 

10 17. A system as claimed in any of claims 11 to 16, wherein 
the relaxionship between che input and output 
frequencies or frequency ranges is given by the 
following 

15 where t, denotes the range of frequencies of the output, and 
4^ and z,^^ denote rhe ranges of frequencies produced by the 
Nth-order and { N -1) th-order nonlinearities, and 



when 
when 



nb 



(a+b) 

nb 
(ai-b]' 



na 

(a + b) J 
na 



<1 



^1 



L (a-^b) J 
n = N and N - 1 
where* [. ] relates to or means take the integer part, 



20 



■ ^ na 
"L {a + b) 



+ 1 



I,, -[na~k{a + b) ,nb-k(a+-b)] for k = 0,...,i' - 1/ 
I, =[0,nb-i*(a + b)], 



and the frequencies of the signal to be processed are in the 
25 range [a>b] and given [a,b] and the required output 
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frequency range fy^ the system further comprises the means 

for derermining rhe smallest: N frora the relationship above 
for the said generalised non- linear system which can bring 
about the specified frequency domain energy transformation. 

5 

13. A aystem as claimed in claim 17, wherein, with N having 
been determined and Kj,,n= 1,,..,N , being given a priori, 
the system further comprises the means: 

for Tzaking N as N and for decsrmixiing the coefficients 
10 of the rime or spatial domain model of rhe generalised 

non-linear system from the description for the system 
output spectrum given in terms of the spectrum of the 
input signal and the coefficiencs of the time or spatial 
domain model of the generalised non-linear sysrem. 

15 

19, ^ A system as claimed in claim 18, further comprising 
means for substituting j w^,"-/ jw^) given m (C13) 

into (C12), and substituting the resultant expression 
for YJjw) into {Cll) to obtain the description for rhe 
20 system output spectrum j.n terms of a function of the 

spectrum of the :i-npat signal and rhe coefficianrs of 
the time or spatial domain model of the generalised 
nonlinear system, 

25 2 0 . A computer progoram product for designing a non-linear 
system for transferring energy from a time or spatial 
domain input signal having a first spectrum at a first 
pre-dererminable frequency or range of frequencies to a 
time or spatial domain output signal having a second 

30 spectrum at a second pre-determinable frequency or 

range of frequencies, the said product comprising, a 
computer readable storage medium comprising: 
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computer program code means for identifying or 
specifying the first spectrum of a zime or spatial 
domain inpur signal from which energy is to be 
transferred/ 

5 coKiputer program code means for specifying the second 

spactriira of a time or spatial domain output signal to 
vhich said energy is to be "transferred, and 

compu-cer program code means for calculating, using a 
frequency domain description of the output sxgnal/ for 

10 example, the output specmuu, expressed in term^ of a 

frequency domain description of the inpur and 
- - coefficients of a time or sparial domain description of 
a generalised non-linear system^ zlne coefficients of a 
time or spatial domain description of said generalised 

15 non-linear system m order to give effec" to the energy 

t3;ansf er . 



21. A computer program product as claimed m claim 20, 
further comprising 

20 computer program code means for select:ing a time or 

spatial domain description of the generalised non- 
linear system; 

computer program code means for determining or defining 
a frequency domain description of the cime or spatial 
25 domain input for the generalised non- linear system; and 

compter program code means fax determining or defining 
the frequency domain description of the output of the 
generalised non-linear system expressed in terms of the 
frequency domain description of said input signal and 
30 the coefficients of the time or spatial domain 

description of a generalised non-linear system. 
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22. A compucer program producr as claimed in either of 
claims 20 or 21^ wherein the frequency domain 
description of the input signal is U(jw), the time or 
sparial domain description of the generalised non- 
linear system is given by the generalised NAKX model 



n-i 



where 



Y.i^) = t t l..,)lly(k-ljriu()«-lj (C181 

with 

Kn t<r 

10 p + q = n, 1,=1,...,K„, i = l,,.„p + q, and ^ = Z'" S 



the frequency domain description of the output of the 
generalised non-linear system is given by 

Y(jw)=X^nt3W) ^^^^^ 

15 where 

N is the maximum order of dominant system nonlinearities^ 



20 denotes an integration over the nth-dur.ensional hyper-piane 
Wi+/-->+w, = w , and jw,,--v3W^) , n = lM^./N are generali^ftd 
frequency response functions of the non- linear system* 

23. A computer program product as claimed in any of claims 
25 20 to 22 y further comprising 
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computer program code means for determining a mapping 
between the i:inie or spatial domain description of the 
generalised nonlinear system and the frequency domain 
description of the generalised nonlinear system* 



24. A computer program product aa claimed in claim 23/ 
wherein the mapping from the time or spatial domaz,Ji 
description of the generalised non-linear system to rhe 
frequency domain description of i:he system is given as 



{C21) 

wtiere 

H,^ljw,,--v]wJ= jH,(jw,,-v3WJH„,,,^,(jw,,,,-,3w^)exp[-j (w, +,-.-,-fw J ij 

(C22) 



25. A computer program product as claimed in any of claims 
2Xi to 24, further comprising 

computer program code means for defining or determining 
a general relationship between the input and output 
20 frequency or frequency ranges of the generalised non- 

linear system. 



26, A computer program product as claimed in any of claims 
20 to 25, wherein the relationship between the input 
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and output frequencies or frequency rar.ges is given by 
the following 

where denotes the range of frequencies c£ the output:, a.nd 
t,^ and t,^^ denotes the ranges of frequencies produced by 

the Nth-ord«r and ( N -1) th-order nonlinearinieS; and 



when 



(a + b)' 



na 

(a^b) 
na 



<1 



n=N and N-1 
where [.] relates no or laeans take the integer part, 



+ 1 



10 =[na-k(a + b) ,nb-'}c(a + b)] for k = 0,..,,!' - 1, 

X^. ={Q,nb-i'(a4.b)], 



(C24) 



and rhe frequencies of the signal to be processed are m the 
range [a/b], given [a^b] and the required output frequency 
15 range t^, the computer program product further comprises 

computer program cod© itiearis for detennining the smallest N 
from the relationship above for the sa^id generalised non- 
linear system which can bring about the specified frequency 
domain energy transformation, 

20 

27. A computer program product as claimed in claim 26, 
wherein, with N having been determined and K^,n = l,,**,N , 
being given a priori, said produce further comprises 
computer program code means far taking N as N and 
25 determining the coefficients of the time or spatial 

domaxn model of the generalised non-linear system from 
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takiin-g N as N and decermir,ing the coefficients of the 
time or spatial domain luodel cf che ganeraiised non- 
linear systeni ^zom tha descripT^ioa far the system 
outpu-t apftctrum ^ivsn Lzi tairias of th^ spactrum of ^ha 
input: signal atr.d the coef f-ciaans Qf th6 tisaa or 
spatial doiuain model of tiae said generalised non-iiaear 
system* 



28. h computer pjrogjr^ product as ci*irtiad iri. any of claiHLS 
:q 20 to ^1, further comprising coo^utar progratr. code 

- means foe substitutxug Hjjw„-/jvJ given in ;C2l) ir.to 
(020^/ ^nd s'obscituting the resultar.t egression for 
Yjjw) into (C19) to obtaiu Lha descriprion ror t]:ie 
nyr^-^iTn output spsctrum In f.^Tins, of a furicricr: r.he 
15 spectsrum o£ the inpuc signal and th^ coefficieuts of 

the tinie or spatial dcaain model ct zhe ganerAliacd 
nonlinear system. 



20 2d. A rion^lineir systera which can transfer energy from a 
time or , spacial domaln^ Input signal having a first 
3pec:trum at a first p«-dat^rmir.ahle frewency or range 
o£ frequencies to a tima or spatial domain output 
signal having a second spectruici at a second pre- 

25 datermlnaTile j:requency or rsui^e a£ frequanoies, said 

syateaa comprising 

means Cor identifying the first spectrum of sha time or 
spatxal d««*in input signal from which energy is to b« 

means for specifying the second spectrum of the tin^e or 
spatial domain output signal to which said energy is to 
be transf err edy and 



30 
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means for giving effect to rhe ervergy transfer usir-g 
corf f icieacs o£ a time or spacial domain description, of 
a gar,eralised noa-Iinear systerci, said coefficients 
liatring been calculatsd using a frstjuency doraain 
44ac5;iptioa q£ 3a.id Au.^put ^igaai, far swsmple, thid 
output spectruis/ expressed in tianna of a frequency 
domain description of said input signal and 
coefficiants of a time or spatial doiiiain description of 
a g-eneraliaed ncn-linaar system. 



30. ' K non.'-Iin^a.r systsnt to tzrans far ring- energy from a tim$ 

or spatial dcrtiain input signal having a first specrrum 
at a first pre-determinable frequency or range of 

15 -'re.q^iencias to -- rrr ap;3r.ial dotiain output signal 

having a sacond spactrum at a second pre-detenciinal^le 
frsquancy or ^ range of frequencies^ said system 
comprising means for iniplemenrinq a cnethod as clainied 
in »ny of claims 1 to 9* 

20 

31, A non-linsar systto for transferring energy froirt a time 
or spatial domain input signal having a first spectrum 
at a first pre-detenninable frequency or range of 
fraqusncies to a time or sparial donia.in output signal 

25 having a second spactrua. at a second pxe^determinabia 

frequency or rstng'e af frequencies, said non^linaar 
system coiKprising a data processing system as claimed 
in arty of claias 11 to 

3G 32. An article of manafacrure coi&prising a computer usabls 
medium with computer readable program code roeana 
embodied in the medium for designing a nan- linear 
system for tr*rts£»r«ing enawgy from a ticiie or spatial 
domain input signal having a first spectrum at a' first 

35 pre-detanainalDie frequency or range of frequencies to a 

S4 
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time or spatial domain output sigr^ai having a second 
spectruitL at a sacand pre-datsrmiriabie frequency or 
range of frequencies, the coiuputer readable prograip. 
code means in said cicticie ggip:pri3ing; 

computer readai^le program code means for identifying 
the first spsctrim of tim« or spatial cioraain input 
signal Srom which energy t;q be tritnafczrredy 

<2otr^ut*r raajlable program code maans far specifying the 
second spectrum of -a time or apatisii domain outpuc 
signal to which said energy is to be transferred, and 

computer * ceadabie program code raeans for c^icalating, 
using a frequency domain description o'f the output 
signal, for example, zhe output spectruni, expressed in 
terms cf: a tL-equency domain descrtf-rxor- the ; :ipv v 

and coefficidnca of * tiiae i^r spatial doir^in 
description of a generalised non-linear sysrem, the 
coefficients of i tiae or spatial domin description oz 
said generalised non-linear systeicv in order co give 
effect: to the energy transfer. 
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NONLINEAR SYSTEM, METHOD OF DESIGN THEREOF AND COMPUTER 

PROGRAM PRODUCT 

ABSTRACT OF THE DISCLOSURE 

The present invention relates to nonlinear systems and methods of design thereof 
in the frequency domain. Typically, conventional linear filter design involves attenuating 
signals at firequencies which are not of interest and dissipating the energy at those 
fi-equencies as, for example, heat or sound. However, in most systems, it is not always 
convenient to design a linear system or design a system solely with energy attenuation in 
mind. Therefore, the present invention provides a nonlinear system and method of design 
thereof in the fi-equency domain which can be used to transfer energy at a first pre- 
determinable frequency, or firequency range, to a second pre-determinabie frequency, or 
fi:equency range. Using the method of the present invention a nonlinear system can be 
developed which can meet given energy transfer requirements or a nonlinear system can 
be designed which can alter the transfer fiinction of an existing nonlinear or linear 
system. 
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p- ^ The power Spectral Density of y ^ 
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Given a sampled series {u(k)} of the signal 
u(t) tp process and the sampling interval T. 



Given a range [c,d] of frequencies specified 
for the output components of the filter to design 



Given a spectrunrr Y*(]w) specified for the spectrum 
of the filter ou^)ut over the frequency band [c,d] 



3002 



3004 



3006 



Calculate U (jw), the spectrum of the input, using 
a fast Fourier Transformation algorithm to yield 
2^t 



..0 M 

2 



where M is the length of the sampled series {u(k)) 



3008 





r ... 


Evaluate the range [a,b] of 
the computed 


the frequencies in u(t) from 
nput spectrum 



3010 



Set 



n^1, Nq=0~|| 



3012 



Perform the toliowing computations 
na 



+1 



Ij=[na-i(a+b). nb-i(a+b)] 

for i^O t*-1 

L^[0,nb-r(a+b)] 



© 

Fig. 30 
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Evaluate the frequency range f y produced 



by me ntn-order nonlinear output as 
ijli when 



i=1 

y Ij . when 
i=0 



(a+b) 
na 

(a+b) 



na 

(a+b) 



na 



<1 



>1 



where ( x ) denotes the integer part of x. 



Whether Nq=0?\ 3018 
Yes ^ No 




3016 



Whether [c,d] has completely failen into L=fy Uf w 

^ 'n n^ 

302r- ■ — 

N=n I 3027 




SetK^=1,for n=Ng N 



© 



3028 



Calculate n 

J exp [-j(w^ +w^lp) n ^ (j^i) 

numerically at w=2^ i/MT i=^0 M/2, 

fori^ =1 ln=^'-''^n' a^^d n-N^^.^N 
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Compute 
round 


"cMT 


and \^ ^(omd 


dMT] 


where round 


(x) means to take the 




integer nearest to x 




\ 
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Based on equations 



f exp [-j(w^ 1 ^ + +Wf, 1^ ) jQ U Gw^ ) d o 

.... PIT . i=1 



Compute 

using at least square routine, 



w 



3034 



Design a conventional linear filter H(q'*'^) to improve the 
above design result within the frequency band [c,d] and 
attenuate the responses over the frequency bands which 
are beyond [c,d]. 



3036 



set Kpj 

forn=Np,...,N 



Check the frequency domain r^ponse 
of the filter ^ 

over the frequency band [c,d] 



3038 



^ No 




Whether the response, in comparision 
with the specified response Y* (jw), 
is okay? 

Yes 



The design is finished and the discrete 
version of the desigried nonlinear system is 

y(k)-H{q-1)=y y Co,(l^ in),IIu(k-ip 



n=H 



'0 iv'n 



=1 
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Given a sampled series (u(k)} of the signal 
u(t) tp process and the sampling intervai I 



Given a range [c,d] of output frequencies 



Given a bound Y^V) specified for the spectrum 
of the filter output over the frequency band [c.d] 



3102 



3104 



3106 



Caicutate U (jw), the spectrum of the input, using 
a fast Fourier Transformation algorithm to yield 

where M is the length of the sampied series {u(k)} 



3108 



Evaluate the range [a,b] of the frequencies in u(t) from 
the computed input spectrum 



3110 



Setn=:1.NQ^0 [ g^^j 



Perform the following computations 
na 



(a+b) 



+1 



Ij=[na-i{a-i-b), nb-i(a'fb)] 

fori^O i'-t 

L-[0,nb-r{a+b)] 



■0 

Fig. 31 
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Evaluate the frequency range fy produced 
Dy the ntn-ofder nonlinear output as 



i=1 

,Uli -w^- 



(a+b) 
na 



(a+b) 
na 



<1 



>1 



(a+b) 

where [ x ] denotes the integer part ot x. 
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Whether Nq=0?J 3118 




Calculate | Uh../|U (iw)| , the nth-order convolution 

integration on |U (jw)[ , numerically to yield 
| U^^,/|U {i2n i/MT )| , i-0 M/2, for n= 
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Compute 






'dMt 




round 


"clvlT 


and =round 










. 271 J 




wnere round 


(x) nneans to take the 




integer nearest to x 





3132 



Based on equations 



U 



m 



Compute ^ ^ ° 

Cn. n=No N. 

using at least squares routine under the constraint that 
the results must be positive 



3134 



Select NARX model parametres 



for DssMq N, under the constraints 



n=NQ N 



and other constraints specified by specific 
design problems 



3136 



Check the frequency donaain response 
of the filter ^ K„ ^ 

y(k)=H(q-1)=£ x^i^ 



n=No I 1^=1 
over the frequency bands beyond [c,d], 
which Should, ideally, be negligible 



3138 




Whether the frequency response 
beyond the band [c,d] is ok? 
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(f) 



The design is finished and the discrete 
version of the designed nonlinear system is 
N '^n n 



3142 




Design a conventional linear band pass filter 
H ( q"'') to attenuate the responses over tfie 
frequency iaands which are beyond [c,d] 



3144 



Complete tne design by constructing a nonlinear 
filter the discrete version of which is 

N "^n n 
y(k)=H(q-l)=2^ ^ CQ,{li..,..ln),ri'J(k'l.) 
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p^ase Angle of Output before Linear Phase Rltering 
over Considered Frequency Band 
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Phase Angie uf Applied Linear Phase Filler 
over Considered Frequency Band 




Frequency (rads/s) 
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Attorney's Docket No.: 12304-002001 
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COMBINED DECLARATION AND POWER OF ATTORNEY 

As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated below next to my name. 

I believe I am the original, first and sole inventor (if only one name is listed below) or an original, first and joint 
inventor (if plural names are hsted below) of the subject matter which is claimed and for which a patent is sought on the 
invention entitled NONLINEAR SYSTEMS. METHOD OF DESIGN THEREOF AND COMPUTER PROGRAM 
PRODUCT , the specification of which: 

[] is attached hereto. 

[X] was filed on September L 2000 as Application Serial No. 09/623,281 and was amended on 

[X was descnbed and claimed in PCT International Application No. PCT/GB99/0Q550 filed on 2 March 
1999 . 

I hereby state that I have reviewed and understand the contents of the above-identified specification, including 
the claims, as amended by any amendment referred to above. 

I acknowledge the duty to disclose ail information I know to be material to patentability in accordance with 
Title 37, Code of Federal Regulations, §1.56. 

I hereby claim foreign priorily benefits under Title 35, United States Code, §1 19 of any foreign application(s) 
for patent or inventor^s certificate or of any PCT international appUcation(s) designating at least one country other than 
the Umted States of America listed below and have also identified below any foreign application for patent or inventor's 
certificate or any PCT international application(s) designating at least one countiy other than the United States of 
America filed by me on the same subject matter having a filing date before that of the apphcation(s) oi which priority is 
claimed: 

Country Application No. Filing Date Priority Claimed 

WIPO PCT/GB99/00550 2 March 1999 [x] Yes [] No 

GreatBritam 9804412.6 3 March 1998 M Yes [] No 

I hereby appomt the following attorneys and/or agents to prosecute this application and to transact all business 
in the Patent and Trademark Office connected therewith: 

Eric L. Prahl, Reg. No.^^^ Frank R. Occhiuti, Reg. No, 35,306 

John F. Hayden, Reg. No. 1LM0_ Y. Rocky Tsao, Reg. No. MJ1S2 

Address all telephone calls to ERIC L. PRAHL at telephone number (617) 542-5070. 

Address all correspondence to ERIC L. PRAHL at: 

FISH & RICHARDSON P.C. 
225 Franklin Street 
__B QSton, MA 021 10-2804 

I hereby declare that all statements made herein of my own knowledge are true and that all statements made on 
mformation and belief are believed to be true; and farther that these statements were made with the knowledge that 
willful false statements and the like so made are punishable by fine or imprisonment or both, under Section 1001 of 
Title 18 of the United States Code and that such willful false statements may jeopai'dize the validity of the application or 
any patents issued thereon. 
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Full Name of Inventor: 

Inventor's Signatui'e; 
Residence Address: 
Citizenship: 
Post Office Address; 
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Great Britain C^ff}/^ 
Great Britain 
16 Grove Road 
Totley, Sheffield S17 4DJ 
Great Britain 



Date: 



Full Name of Inventor: ZI QIANG LANG 



Inventor's Signature: 
Residence Address: 
Citizenship: 
Post Office Address: 



China 
105 Hawley Street 
^^Sh^eMjSl ^A 
Great Britain 
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